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ABSTRACT 

We analyze the physical conditions in the low-ionization component of starburst outflows (in con- 
trast to the high-ionization wind fluid observed in X-rays), based on new Keck/LRIS spectroscopy 
of partially resolved absorption troughs in near-ultraviolet and optical spectra of Ultraluminous In- 
frared Galaxies. The large velocity width and blueshift present in seven, atomic transitions indicate a 
macroscopic velocity gradient within the outflowing gas. The separation of the Mgn 2796, 2803 (and 
Fen 2587, 2600) doublet lines in these data constrains the gas kinematics better than previous studies 
of the heavily blended Nai 5892, 98 doublet. The identical shape of the Mgn 2796 absorption troughs 
to that of the normally weaker transition at 2803A (after accounting for emission filling) requires both 
transitions be optically thick at all outflow velocities. The fraction of the galactic continuum covered 
by the outflow at each velocity therefore dictates the shape of these absorption troughs. We suggest 
that the velocity offset of the deepest part of the troughs, where the covering factor of low-ionization 
gas is near unity, reflects the speed of a shell of swept-up, interstellar gas at the time of blowout. In 
a spherical outflow, we show that the fragments of this shell, or any clouds that expand adiabatically 
in rough pressure equilibrium with the hot wind, expand slowly relative to the geometrical dilution; 
and the covering fraction of low-ionization gas decreases with increasing radius. Our measurement of 
a covering factor that decreases with increasing velocity can therefore be interpreted as evidence that 
the low-ionization outflow is accelerating, i.e. absorption at higher velocity comes from gas at larger 
radii. We also present measurements of Cf(v) in 4 species, place an upper limit of n e < 3000 cm~ 3 on 
the density of the outflowing gas, and discuss lower limits on the mass outflow rate. 
Subject headings: galaxies: starburst — hydrodynamics — infrared: galaxies — intergalactic medium 
— ISM: evolution — line: profiles 



1. INTRODUCTION 

Galactic winds are a key ingredient in cosmological 
models of galaxy evolution. They transport the nucle- 
osynthetic products of stars into galaxy halos and the in- 
tergalactic medium (IGM), thereby shaping enrichment 
and effective yields. The amount of material ejected from 
low mass galaxies can significantly reduce the gas mass 
available for star formation, resulting in lower global effi- 
ciencies for star formation that flatten the faint end of the 
galaxy luminosity function, relative to that of the halo 
mass function (Somerville & Primack 1999). Some heat- 
ing mechanism is apparently required in very massive ha- 
los in order to make the number of galaxies cutoff more 
steeply with luminosity than the halo mass distribution 
does with increasing halo mass. Simulations generally 
invoke some type of feedback from accreting supermas- 
sive black holes, although only the energy scale required 
distinguishes this feedback from supernova-driven, galac- 
tic winds. With mass-loading factors that scale in- 
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versely with galaxy mass, feedback prescriptions, star- 
burst winds account for the metal enrichment of the 
intergalactic medium (IGM;e.g., Oppenheimer & Dave 
2006). 

The scaling relations used for feedback in cosmolog- 
ical simulations are largely motivated by empirical re- 
sults. Nearby starburst galaxies present hot, extraplanar 
gas, heated by supernovae (Dahlem etal. 1998; Martin 
1999; Strickland et al. 2002; Martin, Kobulnicky, & Heck- 
man 2002). Outflow speeds for low-ionization gas have 
been measured for much larger samples of nearby galax- 
ies (Heckman etal. 2000; Rupke, Veilleux, & Sanders 
2002; Martin 2005; Rupke, Veilleux, & Sanders 2005b; 
Sato et al. 2009) and composite spectra of distant galax- 
ies (Shapley etal. 2003; Weiner etal. 2009). In the stan- 
dard dynamical model for the acceleration of this gas, the 
kinematics of the low-ionization gas reflect properties of 
a hotter, energetically dominant wind. The hot wind 
is not a necessary component of the outflow when radi- 
ation pressure (on dust grains) accelerates the outflow. 
Outflow components at vastly different temperatures can 
only be observed from the X-ray to the infrared for rela- 
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tively nearby systems. At redshifts z>0.15, observations 
of winds typically detect low-ionization absorption lines 
and, less commonly, cover strong transitions of Ovi , NV, 
SilV, and CIV. Determining what physical information 
can be reliably extracted from the absorption lines is of 
paramount importance for understanding any evolution 
in galactic wind properties over cosmic time. 

For nearby starbursts, the Nai 5892, 5898 lines are the 
most commonly detected lines from outflows owing to 
the aperture advantage of ground-based telescopes over 
satellites. The majority of luminous starburst galax- 
ies present a strong, blueshifted Nai absorption trough 
(Heckman etal. 2000; Rupke etal. 2002; Martin 2005, 
2006; Rupke etal. 2005b; Sato etal. 2009); but the de- 
tected fraction in Nai drops to about half for the nearby 
dwarf, starburst galaxies (Schwartz & Martin 2004). Al- 
though the stronger line, Nai 5892, is usually saturated 
based on its strength relative to the weaker line at 5898A, 
the absorption troughs are rarely black. These previ- 
ous studies have concluded that the Nai absorbers do 
not completely cover the optical continuum source. The 
blended absorption troughs have been fitted with the 
minimum number of velocity components, typically from 
1 to 3, required to describe the velocity asymmetry (e.g. 
Martin 2005, 2006; Rupke etal. 2002, 2005a); and lit- 
tle has been written about the variation of physical pa- 
rameters with velocity in the outflow. In contrast, the 
very broad absorption troughs identified as outflows in 
active galactic nuclei (AGNs) spectra require a velocity- 
dependent covering factor (Arav et al. 1999a, 2001a; de 
Kool et al. 2002; Gabel et al. 2003; Scott et al. 2004). Re- 
solving absorption troughs in starburst spectra, and com- 
paring different transitions, can lead to new insight about 
the nature of the low-ionization component of outflows. 

In this paper, we present Keck/LRIS spectroscopy 
of redshift z ~ 0.25 galaxies that partially resolves 
the Nai 5892, 5897, Mgi 2853, Mgn 2796, 2803, and 
Fell 2587, 2600 absorption troughs. We choose Ultra- 
luminous Infrared Galaxies (ULIRGs) due to their high 
Nai outflow fraction and selected systems at redshifts 
that provided Nai and near-ultraviolet spectral cover- 
age. The targets are among the most luminous star- 
bursts in the local universe and are classified as ULIRGs, 
log Lir > 12. Their activity triggers by gas inflow 
induced by a recent, or on-going, merger (Borne etal. 
2000). The wider separation of the Mgn doublet lines, 
768 km s _1 , relative to the Nai doublet, allows direct 
comparison of the A2796 and A2803 troughs, hereafter 
the blue and red troughs respectively. At any location 
in the outflow, the blue transition has twice the opti- 
cal depth of the red one; and the covering fraction will 
necessarily be identical. In contrast to the a-element 
enrichment measured in the hot wind (Martin, Kobul- 
nicky, & Heckman 2002), the relative abundance of Fe 
to Mg could be as high as solar in the low-ionization 
gas, which may be primarily entrained interstellar gas. 
The optical depth in the Fen line would then be just 1.8 
times lower than that in Mgn 2803. However, the much 
weaker Fen 2587 line, r (2600) = 3.5r (2587), might be 
optically thin when the other transitions are saturated. 
Whether Fell 2587 or Nai 5898 proves to have lower op- 
tical depth, and therefore provides the better measure- 
ment of ionic column density, depends on the relative 
ionization corrections for Fen and Nai . 



Our presentation is organized as follows. Section [2] de- 
scribes the new observations, derives emission-line red- 
shifts, and estimates the sensitivity to outflowing gas. 
On the first reading, we recommend skipping directly to 
§ [H where the absorption trough measurements are first 
presented in a model-independent fashion. Section 13.21 
then uses the independent measurements of the trough 
intensity Ib(v) and Ir(v) to solve directly for optical 
depth tr(v) and covering fraction C/(v) as a function of 
outflow velocity. Since this approach cannot be applied 
to the majority of the absorption troughs, we introduce a 
Xy-fitting method in § 13.31 It differs from previous anal- 
yses in terms of the set of physically motivated priors 
adopted. The advantages are quantitative measurements 
of the velocity dependence of the covering factor in four 
species and direct comparison of the limits on gas col- 
umn density from different transitions. In this way, we 
establish the relation of the outflow component probed 
by near-UV resonance lines to that seen in Nai . The 
connection enables comparison of wind properties over a 
large redshift range. Our discussion in Section 0] focuses 
on interpreting the velocity-dependence of the covering 
fraction. The main results and their implications are 
summarized in Section [5] 

Oscillator strengths and vacuum wavelengths are taken 
from Morton (1991, 2003) for the optical and near- 
ultraviolet transitions. We adopt a cosmological model 
with H = 70 km s" 1 Mpc" 1 , Q = 0.3, and ft A = 0.7 
throughout the paper. 

2. KECK OBSERVATIONS OF OUTFLOWS IN THE 
NEAR-ULTRAVIOLET AND OPTICAL 

We obtained (rest-frame) near-UV and optical spec- 
troscopy of redshift ~ 0.25 starburst galaxies with LRIS 
(Oke etal. 1995; McCarthy etal. 1998) on Keck I. The 
blue channel of LRIS is one of the few optical spectro- 
graphs that is efficient down to the atmospheric cut-off 
in the blue, allowing observations of Mgll at z> 0.15 
(for targets directly overhead). The dual-beam design 
offers simultaneous coverage of Nai absorption and sev- 
eral emission lines in the red channel. Due to the low 
density of bright galaxies at these redshifts, the spectra 
must be obtained one at time. 

Galaxies were chosen from the IRAS 1 Jy survey 
of ULIRGs (Kim & Sanders 1998), which contains 118 
ULIRGs with Fqq > 1 Jy, based on apparent magnitude 
and airmass at the time of our runs. Ultraluminous In- 
frared Galaxies are extremely rare at the present epoch 
but offer the closest local analog of the IR-luminous 
galaxy population that dominates the cosmic star forma- 
tion rate density at redshifts greater than 0.7 (Le Floc'h 
etal. 2005). The high dust content absorbs much of the 
stellar luminosity and emits this energy as thermal radi- 
ation in the infrared. 

Table [1] lists some properties of the galaxies whose 
spectra are presented in this paper. The infrared lumi- 
nosities of the galaxies observed range from log Ljr = 
12.31 — 12.81 with mean logL]R = 12.58, which is 
more luminous than the mean of the Martin (2005, 
2006) sample. The AGN fraction is higher in brighter 
ULIRGs (Veilleux etal. 1999); and two of the five 
ULIRGs in the sample are spectroscopically classified 
as AGN (2 LINERs and 2 Sey 2's). The more recent 
Kewley et al. (2006) division of these optical-emission- 
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line ratios into excitation types classifies FSC 0039- 
13 as a composite starburst-AGN, FSC 1009+47 as a 
Seyfert 2, FSC 1407+05 at the Seyfert/LINER bound- 
ary, FSC 1630+15 at the LINER/Composite boundary, 
and FSC 2349+24 as a Seyfert 2. Among our five targets, 
FSC 2349+24 clearly presents the strongest AGN signa- 
tures, but AGN contribute less than half of the bolomet- 
ric luminosity in the other four ULIRGs (Veilleux et al. 
2009a). Among our five targets, FSC 2349+24 clearly 
presents the strongest AGN signatures. Rupke et al. 
(2005d) found little difference in Nal wind kinematics be- 
tween HII, LINER, and Sey 2 ULIRGs; only the Sey l's 
showed faster outflows. Hence, we expect these outflows 
are driven primarily by the starburst rather than the 
AGN. The far-infrared colors are all cooler than the stan- 
dard demarcation, F(25/xm)/F(60^m) > 0.2, used to 
identify warm ULIRGs (Kim et al. 1998), which also sug- 
gests the total luminosity is starburst dominated. The 
total infrared luminosities correspond to upper limits on 
the star formation rates from SFR ~ 207 — 653 M Q yr _1 
or ~ 352 — 1110 M Q yr _1 , respectively, for the Chabrier 
or Salpeter initial mass functions. 

2.1. Data Acquisition and Reduction 

Longslit spectra were obtained with LRIS (Oke etal. 
1995) and LRISb (McCarthy et al. 1998) on 2004 January 
26, 2004 March 16-17, 2007 October 6, and 2007 Novem- 
ber 1. Clouds and high humidity limited the exposure 
times on all of these nights. We present the five highest 
quality spectra obtained from the combined data. 

The position angle of the longslit was selected to cover 
both nuclei of FSC 1009+47 and FSC 2349+24. The at- 
mospheric dispersion compensator, which was installed 
on LRIS in spring 2007, was used for the 2007 obser- 
vations. Care was taken with the 2004 observations to 
observe the targets when the slit PA was near the paral- 
lactic angle. 

The dual beam spectrograph was configured with the 
D560 dichroic, a 1200-1 grism blazed at 3400A in the blue 
arm (LRISb) and a 1200-1 grating blazed at 7500A in the 
red arm. The grating tilt for each target was tuned to 
cover the Ha + [Nil] emission lines as well as the Nal 
absorption and Hel 5876 emission line. The blue spectra 
cover the Mgn 2796, 2803 doublet, the Mgi 2853 line, 
and the Fell 2587, 2600 doublet . The slit width was 
chosen based on the atmospheric seeing. The resolution 
(for a source filling the slit) ranged from 110 km s _1 on 
the 2004 March run to 160 km s^ 1 on the 2007 November 
run and is listed in Table [2] The size (FWHM) of the 
galaxies along the slit exceeded the slit width. 

Fixed-pattern noise was removed using software scripts 
that called IRAF tasks. 2 The blue spectra were flat- 
fielded with twilight sky frames, normalized by the sky 
spectrum. The red spectra were flatfielded using an in- 
ternal (spectrosopic) exposure of a quartz lamp. A dis- 
persion solution was fitted to the vacuum wavelengths 
of the arc lamp lines as a function of detector coordi- 
nate. The root mean square error in the dispersion so- 
lution for the blue (red) spectra was 0.05 (0.07) A. Ap- 
plication of a small, additive shift, up to a couple tenths 

2 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 



See Figure 1 on p. 22 

Fig. 1. — Normalized intensity vs. velocity, where the velocity is 
relative to the systemic velocity determined from emission lines. 

of an angstrom, registered the wavelengths of night sky 
emission lines with their values in a telluric-line spec- 
trum (Hanuschik 2003) smoothed to our spectral resolu- 
tion and transformed to vacuum wavelengths (using the 
Edlcn formula). We attribute these corrections to shifts 
in the dispersion solution with airmass and rotator angle. 
The correction to the Local Standard of Rest (LSR) was 
computed for each observation using the IRAF task RV- 
CORRECT. All the offsets were less than 30 km s" 1 . 
Since we are only concerned with relative velocites, the 
corrections to LSR were not applied to the data. 

We rectified the two-dimensional spectral images, us- 
ing the dispersion solution and traces of a standard star 
stepped along the longslit, and then extracted an in- 
tegrated galaxy spectrum for each target. These spec- 
tra have SNR ^5—10 per pixel as shown in Ta- 
ble [21 Additional, lower quality spectra were extracted 
for the distinct continuum sources within FSC 1009+47, 
FSC 1407+05, and FSC 2349+24. After substraction of 
the median sky intensity at each wavelength, significant 
residuals from strong night-sky emission lines remained 
in the red spectra near the Nal line in FSC 1009+47, 
1407+05, and 1630+15. Variance spectra were extracted 
for each target prior to sky subtraction and flux calibra- 
tion. The variance vectors were later scaled by a multi- 
plicative factor that made the uncertainties in the inten- 
sity consistent with the measured standard deviation in 
the extracted, target spectra. 

Observations of multiple, standard stars determined 
the relative sensitivity with wavelength. The data were 
flux calibrated using this sensitivity function and then 
normalized by a fitted continuum. The error in con- 
tinuum placement is ncligiblc around Nal , Mgi , and 
Mgn ; but the blending of stellar absorption lines washes 
out the true continuum level at shorter wavelengths. To 
identify bandpasses near Fen 2587, 2600 that likely reach 
the true continuum level, we compared high-resolution, 
synthesized spectra models of stellar populations to 
copies smoothed to 100 km s _1 resolution. We fitted 
a first or second order cubic-spline through these band- 
passes as well as the broad bandpasses near the Mgi and 
Mgn lines. The resulting error in the continuum level 
near Fen depends on the star formation history. The 
severe blanketing in older bursts requires actual fitting 
of reddened, model spectra. This level of sophistication 
was unneccessary, however, because these near-UV spec- 
tra are bluer than the t > 100 Myr burst models; and 
we confidently rule out old, burst models for the con- 
tinuum. The population synthesis models indicate ei- 
ther continuous star formation or a burst within the past 
100 Myr. Repeated fitting trials indicate the uncertainty 
in the continuum level around Fell was 1 to 5%. 

For a typical line width of 470 km s _1 FWHM, 
we detect a rest-frame equivalent width W r (5a) w 
1.52A(5/ SNR pzx ) (Av/470 km s -1 ) 5 in the red spec- 
tra at the 5<7 significance level. A typical 5cr 
sensitivity limit for the blue spectra is W r (5a) » 
0.92A{5/ SNR pix ) (Av/470 km s" 1 ) - 5 . 

2.2. Sensitivity of Spectra to Physical Properties 
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We measured redshifts from recombination lines of H 
and He and forbidden lines of singly-ionized N and S. 
These lines all fall in our red spectra. The dispersion 
solution ties both the blue and the red spectra to the 
vacuum wavelengths of night sky emission lines. These 
redshifts agree with the previous measurements of Kim & 
Sanders (1998). Our independent check is important be- 
cause that work used the Nai absorption line in the red- 
shift estimate; and we will show that the Nai kinematics 
differ significantly from that of the recombination lines. 

The sensitivity of the seven metal lines to outflowing 
gas depends on the relative abundances of the elements, 
the oscillator strengths of the transitions, the dust deple- 
tion, and the ionization corrections. In the limit of com- 
parable ionization corrections for all species, the optical 
depth in the Mgi 2853 line would be the largest for solar- 
abundance ratios. The Mgn lines have lower oscillator 
strengths, and the cosmic abudances of Fe and Na are 
lower. Our typical detection limit for Mgi corresponds 
to a minimum hydrogen column density, 

N(HI) = im x 10 17 cm~ 2 (N(M g) /N(M gl)) x 

1 OK 

(5/SNR) ——(Aw/470 km s" 1 ) - 5 . (1) 

1 ~\~ z 

Under the same conditions, the Nai 5898 line optical 
depth would be slightly lower than the weakest Fen line, 
but whether either transition is optically-thin depends 
on the relative ionization correction, where 

tq (Fen 2587) = v(Fell ) 

r (Nai 5898) ' x(Nai )' 1 ' 

The Nai 5898 line can, in principle, provide the best mea- 
surements of an ionic column density in the outflow, but 
interpretation is complicated by the blending with the 
Nai 5892 transition, the ionization correction, and de- 
pletion of Na by dust grains. The Fen transitions with 
the lowest oscillator strengths lie blueward of the atmo- 
spheric cut-off for nearby galaxies, and their measure- 
ment will likely lead to the best estimates of mass column 
density. 

3. MEASURED PROPERTIES OF LOW-IONIZATION 
OUTFLOWS 

Figure [1] shows the absorption troughs on a velocity 
scale. The high spectral resolution, relative to the broad 
troughs, allows us to compare the shape of the absorption 
troughs among seven transitions. The centroids of the 
Ha and [Nil] emission lines lie at the systemic velocity 
by construction. The spectra of all five galaxies present 
broad absorption troughs in the Nai and Mgn dou- 
blets. The spectrum of FSC 0039-13 has the best SNR 
near the Mgi 2853, Fen 2600, and Fen 2587 lines; but 
these transitions are detected in all the spectra except 
for FSC 2349+24. The Fen 2587, 2600 lines lie short- 
ward of A3200 for FSC 2349+24, and the SNR is severely 
compromised by atmospheric attenuation. In Figure [T] 
the two Seyfert 2 galaxies in the sample, FSC 1407+05 
and FSC 2349+24, present Mgn in emission. Some He I 
emission is present near the Nai doublet in all the spec- 
tra and is most prominent in FSC 0039-13. 

Spectra of A, F, and G stars present low-ionization, 
metal lines; but any photospheric contribution to the to- 
tal equivalent width in our spectra is small. The velocity 



offset of the absorption troughs require a non-stellar ab- 
sorption component; and the Mgn troughs for four of 
the galaxies show little or no absorption at the systemic 
velocity. Only the FSC 0039-13 spectrum presents ab- 
sorption at v = 0. We attribute it to the interstellar 
medium in the host galaxy, likely a result of high incli- 
nation relative to our sightline. 

A dominant stellar origin for the resonance absorp- 
tion appears unlikely in these ULIRGs for several rea- 
sons. First, Figure [2] shows Mgn and Mgi equivalent 
widths measured from synthetic spectra, computed with 
the STARS2002 stellar population synthesis code (Stern- 
berg 1998; Thornley et al. 2000; Sternberg, Pauldrach, & 
Hoffmann 2003; Davies etal. 2007) and the UVBLUE 
library of high-resolution stellar spectra (Rodriguez- 
Merino etal. 2005), for a broad range of star for- 
mation histories. 3 Only post-starburst populations, 
with continua dominated by A, F, or G stars, show 
Mgll 2796 equivalent widths of more than a cou- 
ple Angstrom. The large Mgn equivalent widths of 
FSC 0039-13, FSC 1407+05, and FSC 1630+15 are in- 
consistent with continuous star formation; and the other 
two ULIRGS do not lie near the stellar locus. Second, 
the ULIRG spectra lack detectable absorption from ex- 
cited, electronic states of Fell . lines that are prominent 
in the spectra of older stellar populations in Figure [H 
In addition, the synthesized spectrum of the stellar pop- 
ulation shows Mgn lines much broader than those from 
Mgi , although smearing by the motions of stars in the 
galaxy could hide this difference. Previous ULIRG out- 
flow studies estimated the stellar contamination in the 
Nai absorption trough from the equivalent width of the 
excited Mgi triplet at ~ 5200 A (not covered by our spec- 
tra) and typically found a negligible stellar contribution 
(Martin 2005). 

Extracting quantitative information about the absorb- 
ing material in the outflow from the absorption troughs 
requires a physical model. The spectra presented in 
this paper resolve components with intrinsic velocity 
widths of 100 km s _1 or more; and neither thermal nor 
turbulent motion easily explains the extremely broad, 
~ 800 km s _1 , velocity width of the absorption troughs. 
Macroscopic variations in velocity are required and could 
be associated with discrete supershells and their frag- 
ments (Fujita etal. 2009), interstellar clouds over-run 
by the superbubble shock (Cooper etal. 2008), eddies 
formed at the wind - disk interface (Heckman et al. 
2000), or the velocity gradient in a smoothly accelerating 
wind (Murray et al. 2005). These ideas motivate descrip- 
tions with discrete velocity components, where compo- 
nents correspond to individual clouds or shell fragments. 

In Section 13.31 we formally fit all the absorption 
troughs with such velocity components. For unblended, 
doublets, a parametric description of the absorption 
troughs requires fewer physical priors, however; and we 
apply this approach in Section ^. 2l to build intuition. The 
key constraint throughout this modeling is the observa- 
tion that many of the transitions must be optically thick, 
not only at the deepest part of the absorption trough but 
at high outflow speed. To illustrate the robustness of this 

3 All models assumed a Kroupa initial mass function from 1 
to 120 Mq, solar metallicity, and an exponentially declining star- 
formation rate. 
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model-indepdendent statement, we directly compare the 
absorption trough shapes among the different transitions 
in Section I5TT1 
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Fig. 2. — Mgl and Mgn absorption equivalent wid ths f or 5 
ULIRGs compared to simple stellar population models (§ [3TTJ and 
individual stellar templates. The equivalent widths increase with 
increasing age, from logr(j/r) = 7 to 10, for both burst (solid 
squares) and continuous (solid triangles) star formation histories. 
In the individual stellar templates, the EW(Mgll ) is set by the 
effective temperature, while the EW(Mgl) varies significantly with 
the gravity parameter. In § 13.11 we argue that neither the Mgl nor 
Mgn absorption equivalent widths are stellar-dominated in any of 
these ULIRGs. 



3.1. Direct Comparison of Absorption Trough Shape in 
Seven Transitions 

When intervening absorption is studied in quasar spec- 
tra, or the Galactic halo studied in absorption against 
stellar spectra, the angular size of the absorbing gas 
clouds exceeds that of the continuum source. Gas be- 
tween the observer and the light source attenuates the 
continuum by an amount proportional to the logarithm 
of the optical depth in any transition. For a doublet, we 
have 

I B (v) = I e- TB ^ (3) 

Ir(v) = I e- Tniv) . (4) 

For interstellar conditions, the relative optical depth 
between electronic transitions in a single ion from the 
ground state, i.e. of resonance lines, is tb/tr = 
(/BAs)/(/flAij). For the Nai and Mgn doublets, the ra- 
tio of the oscillator strengths of the blue line, fs, to 
the red line, fa, is 2; and the wavelengths of the tran- 
sitions are very similar. Substituion in Eqn. [3] indicates 
that the relative depth of the continuum-normalized ab- 
sorption troughs must be Ib(v) — Ir{v). These rela- 
tions hold provided the clouds completely cover the con- 
tinuum source over a velocity range comparable to (or 
greater than) the spectral resolution. In practice, when 
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Fig. 3.— Unreddened STARS2002/UVBLUE models for a 1 Gyr 
old stellar population of solar metallicity vs. the FSC 0039-13 spec- 
trum (in black) around Fen . The data show very strong Fell lines. 
The weaker Mnll 2576.877, 2594.499, and 2606.462 lines are ex- 
pected to have equivalent width ratios of 1.8:1.4:1.0, respectively. 
Since the upper limit on the Mnll 2606 strength is low, the ob- 
served Fell 2600 trough cannot have a significant contribution from 
an Mnll 2594 line. The spectral quality allows only upper limits to 
be placed on the strength of the excited Fen * lines. Excited Fell * 
lines are prominent in the burst model and present but weaker than 
our detection limit in the continuous star formation model. Mod- 
els were normalized to the data at 3000A. Modest reddening of the 
continuous star formation model, e.g. Ay = 0.5, provides a reason- 
able estimate of the continuum level in the FSC 0039-13 spectrum. 
Burst models with ages greater than 100 Myr are intrinsically too 
red to describe the observed continuum slope. 

transitions with different / values saturate, both troughs 
appear black in noisy spectra. 

Observations of Nai absorption troughs in ULIRGs 
cannot directly measure the relative intensities Ib(v) to 
Ir(v) due to the blending of the two absorption troughs. 
Previous modeling of the combined trough strongly sug- 
gests, however, that the assumption of complete contin- 
uum coverage is not met (Martin 2005,2006; Rupke etal. 
2005a,b). 

3.1.1. Description o/Mgll Absorption Troughs 

The Mgn absorption troughs delineate the outflow 
kinematics most cleanly. In FSC0039-13, the majority 
of the trough equivalent width is blucshiftcd, but the 
minimum intensity of the trough lies near the systemic 
velocity. The intensity minima are blueshifted several 
hundred km s _1 in the other Mgn spectra with little or 
no absorption at the systemic velocity. In FSC1630+15 
and FSC2349+24, the minima in the absorption troughs 

lie at v 400 km s" 1 . The FSC1407+05 spectrum 

presents a broad trough from -200 to -500 km s _1 . 
The Mgn 2796 and 2803 absorption troughs blend to- 
gether in the FSC1009+47 spectrum, possibly due solely 
to the lower SNR of this spectrum. The Mgn absorption 
troughs have residual intensities < 10% in FSC0039-12 
and FSC 1630+ 15 at minimum intensity. The residual 
intensity at line center is low but not zero in the other 
three spectra. The Mgn absorption troughs are notice- 
ably deeper than the Mgl or Nai absorption troughs. 

We can test the condition represented by Equa- 
tions 3 and 4 by comparing the Mgn A2796 and A2803 
troughs in our spectra of FSC0039-13, FSC1407+05, 
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FSC1630+15, and FSC2349+24, only FSC 1009+47 
presents severe blending of the absorption components 
in Figure [T] Comparison of the Mgn 2796 and 2803 
trough shapes illustrates their nearly identical intensity 
at most velocities. The most significant differences ap- 
pear in the FSC2349+24 and FSC1407+05 spectra where 
Mgn A2796 emission fills in part of the A2803 trough. 
The assumption of unity covering factor implicit to the 
apparent optical depth method clearly fails. As seen 
in projection on the sky against the spatially extended 
galactic continuum, the low-ionization outflow does not 
uniformly cover the source. 




-1000 -500 500 1000 



AV (km/s) 

Fig. 4. — Mgl with scaled Mgn profiles are overlaid. To fa- 
cilitate comparison the depth of the corresponding Mgll absorp- 
tion troughs have been reduced by factors of 2.0, 2.5, 2.0, 2.0, 
and 1.0, respectively. The shape of the Mgll and Mgl absorption 
troughs are indistinguishable aside from the contribution of the 
Mgn emission. 

3.1.2. Relative Shape o/Mgl and Mgn Absorption Troughs 

In FigureH] we scale the depth of each Mgn absorption 
trough to facilitate comparison to the shape of the 
Mgl absorption trough. Based on our inspection of 
the data, we argue that the Mgn and Mgl absorption 
troughs trace gas with the same kinematics. 

FSC0039-13:— We reduced the depth of the 
Mgn absorption trough by a factor of two by mul- 
tiplying I c {\) — by a factor 0.5 and adding the 
result to the normalized continuum, / C (A). The ab- 
sorption troughs in Mgn 2796, 2803 and Mgl 2853 are 
essentially indistinguishable after this scaling. The 



Mgn absorption towards FSC0039-13 is detected at 
velocities ~ 100 km s _1 higher than the maximum 
velocities detected in Mgl . The SNR in the Mgl profile 
blueward of -600 km s _1 is inadequate to definitively 
detect the weak absorption expected based on the shape 
of the Mgn 2796 profile, so the Mgl observation is 
consistent with the presence of neutral Mg up to the 
same outflow velocity detected in Mgn . 

FSC1009+47:— The scaled Mgn 2796 and 2803 absorp- 
tion trough is identical to that of Mgl blueward of - 
100 km s _1 out to -600 km s _1 , where Mgn 2803 blends 
with the red part of 2796. The high- velocity 2803 absorp- 
tion makes the 2796 trough appear to be deeper than the 
2803 trough between -100 km s _1 and +150 km s" 1 . 
The 2796 trough is detected to -780 km s" 1 . The 
Mgl trough is cleanly detected to -450 km s~ x ; but com- 
parison to Mgn 2803 suggests Mgl is plausibly present in 
the outflow to at least -700 km s _1 . After the depth 
of the Mgn profile is reduced by a factor of 2.5, the 
Mgl profile is indistinguishable from the 2796 profile at 
high velocity and the 2803 profile at low velocity. 

FSC1407+05:— The Mgl and Mgn 2796 lines in 
FSC1407+05 are both detected to velocities reaching - 
750 km s^ 1 . The Mgll 2803 profile presents emis- 
sion from to +250 km s^ 1 . The Mgll 2796 emis- 
sion is diluted by the Mgll 2803 absorption at the same 
wavelength. The addition of the 2796 emission to the 
Mgn 2803 absorption trough elevates the latter between 
velocities -750 and -500 km s _1 . At intermediate ve- 
locities less affected by blending, -500 to km s _1 , the 
Mgn and Mgl profiles lie right on top of each other after 
reducing the depth of the Mgn troughs by a factor of 
2.0. 

FSC1630+15:— In FSC1630+15, the shape of the Mgn 
and Mgl profiles are the same from to -620 km s _1 af- 
ter scaling Mgn by a factor of 2.0. Maximum absorp- 
tion is offset to -400 km s" 1 . The Mgll 2796 absorp- 
tion trough is detected to -825 km s _1 . The SNR in 
Mgl near these velocities is not adequate to detect the 
scaled Mgn profile, so the velocity ranges of absorption 
Mgl and Mgn are consistent with being the same. 

FSC2349+24:— The FSC2349+24 spectrum shows 
strong Mgn emission. The 2796 emission fills in the 
bluest portion of the 2803 line profile. The 2796 trough 
reaches -650 km s _1 . We compare the Mgll trough to 
Mgl without any scaling. The SNR is not adequte to 
detect the Mgl trough at outflow speeds greater than - 
500 km s" 1 . 

3.1.3. The Mgl Absorption Trough as a Template for the 
Blended Nal Trough 

Since we detect these ULIRG outflows in neutral Na, 
the spectral energy distribution incident on the outflow 
likely presents a Lyman edge. (See the quantitative dis- 
cussion in §2.3 of Murray etal. 2007). Stellar spectra 
have large intrinsic Lyman edges, so it follows that the 
starburst, rather than AGN, likely determines the Na 
ionization of the outflow. We expect a similar situation 
for Mg, which has an ionization potential just 2.5 eV 
higher than that of Na (7.6 eV vs. 5.1 eV). The patchy 
distribution of dust in ULIRGs, however, leaves open the 
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Fig. 5. — Comparison of Nal and Mgl absorption line troughs in 5 ULIRGs. Left: The blue side of the Nal trough, dominated by the 
Nal 5892 transition, compared to the Mgl trough. Right: The red side of the Nal trough, dominated by the Nal 5898 transition, compared 
to the Mgl trough. 
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possibility that the morphological appearance of the stel- 
lar component in the near-UV differs considerably from 
that at 6000A. Absorption troughs present in these two 
bandpasses need not sample the same region of the out- 
flow. 

Before comparing the near-UV and optical trough 
shapes for the first time, it is useful to consider the rela- 
tive optical depths in each transition should the absorp- 
tion occur in the same clouds. Magnesium is 19 times 
more abundant than Na but only four times more de- 
pleted than Na. 4 Singly ionized magnesium, Mgn , is the 
dominant ion of Mg over a large range of density and tem- 
perature. In contrast, the Murray et al. (2007) ionization 
models predict a larger variation in x(Nai ) from object 
to object, allowing a wide range of possibilities for the 
relative optical depth of Mgi and Nai . For comparable 
Mg and Na ionization fractions, i.e. x(Nal ) w x(Mgn ) 
where x(X+ n ) = N(X +n )/N(X), the Nai 5892 optical 
depth and that of the weaker Nai 5898 line will be lower 
than the Mgi 2853 optical depth. 

The blending of the Nai 5892, 98 lines complicates 
the direct, model-independent comparison. In Figure [5l 
we use the Mgi profile as a template for both the 
Nai 5892 (left column) and 5898 (right column) absorp- 
tion troughs. In three of the ULIRGs, the depth of the 
Nai and Mgi absorption troughs are similar at minimum 
intensity. The two ULIRG spectra with strong Hel 5877 
emission present shallower troughs in Nai relative to 
Mgi , which we explain in part by emission filling of the 
Nai troughs. We take account of the Hel 5877 emission 
filling in Section 13.31 and estimate comparable maximum 
outflow speeds in Nai , Mgi , and Mgn . Surprisingly, 
the modeling in Section 13.31 demonstrates that the kine- 
matics of the Nai absorbing gas are not distinghuishable 
from that in the low-ionization UV transitions. 

FSC0039-13:— In the FSC 0039-13 spectrum, we de- 
tect Nai absorption up to -320 km s _1 . Emission 
filling from He I plausibly explains why the Mgi and 
Mgn absorption troughs extend to higher velocity. The 
Mgi trough is nearly twice as deep as the bluer portion of 
the Nai trough (left panel) and three times deeper than 
the redder portion of the Nai trough (right panel). The 
Hel emission contributes to the shallowness of the blue 
half of the Nai trough (relative to Mgi ), but emission 
filling fails to expalin the low intensity of the red half of 
the Nai trough over 1000 km s" 1 from the Hel line. The 
deepest part of the Mgi profile at +90 km s _1 matches 
the minimum in the Nai trough. At the wavelength 
corresponding to the same Doppler shift in the Nai 
5892 transition, we find another local minimum in the 
blended, absorption trough. We postpone discussion of 
the relative intensities of the Nai 5892 and 5898 absorp- 
tion troughs to Section |3~31 because both transitions con- 
tribute to the absorption at minimum intensity. 

FSC1009+47:— In the FSC1009+47 spectrum, the 
Nai absorption trough extends from roughly the sys- 
temic velocity to about -450 km s _1 . At the coarse level 
of comparison allowed by the poor spectral SNR, the 
shapes of Nai 5892 and Mgi profiles match without any 

4 Depletion quoted for cool interstellar medium of the Galaxy 
(Savage & Sembach 1996). 



scaling. Better SNR at Mgn shows the outflow reaches 
-800 km s" 1 . 

FSC1407+05:— The Nai absorption is detected to - 
650 km s _1 , close to (within one resolution element of) 
the maximum absorption velocity detected in Mgi and 
Mgn . The Nai and Mgi absorption trough shapes are 
very similar, with the latter being slightly deeper. 

FSC1630+15:— In FSC1630+15, the Mgi profile pro- 
vides a good description of the Nai absorption trough 
with no scaling of the Nai 5892 or 5898 lines. Absorption 
is detected in both troughs to -000 km s -1 . The SNR's 
of the Nai and Mgi spectra do not allow detection up to 
the maximum velocity detected in Mgn . 

FSC2349+24:— In FSC2349+24, Hel emission fills part 
the blue-end, v < —250 km s _1 , of the Nai absorption 
trough. The Mgi and Nai troughs have similar intensity 
at other velocities. 

3.1.4. Relative Strength of the Fen and Mgn Absorption 

Troughs 

The oscillator strength of Fen 2600 is about 3.5 times 
higher than that of Fen 2587. The similar shape of the 
two blueshifted Fell troughs requires even the Fen 2587 
transition to be optically thick at the highest veloci- 
ties measured in the outflow. Iron has only a slightly 
lower cosmic abundance than Mg. Using cosmic abun- 
dance ratios, the optical depth of Fen 2600 will be 57% 
that of Mgn 2803 for comparable depletion and ioniza- 
tion corrections. We directly compared the shape of the 
Fen 2600, Fen 2587, Mgn 2796, and Mgn 2803 absorp- 
tion troughs below, omitting the two objects with low 
SNR near Fen . We find similar intensity in all four 
troughs at comparable velocity. Section 13.31 presents 
quantitative, but model dependent, fits. 

FSC0039-13:— From -600 km s" 1 to +200 km s" 1 , 
I B ~ Ir for the Fen 2587, 2600 troughs. The Fell 2587 
trough does not present the dip seen in Fen 2600 from 
-800 to -600 km s" 1 . We do not detect the Mnii 2576 
line, so absorption from the weaker Mnii 2594 transi- 
tion does not produce the blue dip in Fen 2600. We 
also exclude blending with excited, photosphcric lines 
based on the location of such lines in the stellar spec- 
tra shown in Figure [3] We attribute the discrepancies 
in the Fen profile shape over a small velocity range to 
a combination of statistical and systematic errors in the 
data, where continuum fitting uncertainties dominate the 
systematic error in this spectral region. The Fen 2600 
trough is identical to the Mgn trough (to within the er- 
ror bars) over the velocity range from -600 km s" 1 to 
+400 km s _1 . The Fen 2587 trough varies in lock step 
with the Mgn trough from -900 km s _1 to +200 km s _1 . 
We conclude that the Fell and Mgn absorption troughs 
sample gas with similar kinematics. 

FSC1009+47:— The Mgn profile for FSC1009+47 pro- 
vides a good description of the blue wing of the Fen 2600 
trough with no scaling. The Fen 2587 trough is slightly 
shallower. 

FSC1407+05:— The blueshifted wing of the 
FSC1407+05 Fen spectra at v < -200 km s _1 shares 
the shape of the Mgn profile, but the Fen troughs are 
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slightly shallower (amplitude reduced by 1.3 in 2600 
and 1.6 in 2587). Both Fell troughs lie well below the 
Mgll troughs from -100 km s _1 to +100 km s _1 . The 
Mgn 2803 profile presents a weak, but significant, emis- 
sion bump from to +200 km s _1 . No Fell emission 
is detected, so the Fell absorption troughs will look 
broader than the Mgll troughs. The Mgn absorption 
trough may share the kinematics of the low-velocity 
Fen trough but be filled in by emission. 

3.2. Parametric Descriptions of the Absorption Troughs 

In the spectral regions where the Mgn 2796, 2803 ab- 
sorption troughs can be directly compared, their inten- 
sities are nearly equal, Ir ~ Ir- Since the galaxy (con- 
tinuum) likely subtends a large angle relative to that 
of an individual shell fragment or cloud, we might ex- 
pect variation in the properties of the absorbing material 
as a function of the spatial and/or velocity coordinates 
in the outflow. A number of methods have been devel- 
oped to parameterize the relative geometry of the out- 
flow and continuum source. We apply the most common 
method to the absorption troughs from FSC0039-13 and 
FSC1630+15, which show no sign of emission filling and 
little blending of the A2796 and A2803 lines. We then 
explain why two other parameterizations fail to describe 
the distribution of outflowing gas. 



F:::ru039-13 




Fig. 6. — Covering fraction vs. Mgll velocity. We binned to 
the spectral resolution and computed Cf(v) from Equations [7] 
For FSC0039-13, the open symbols show the effect of removing 
a symmetric, zero-velocity component. For reference, the solid 
(blue) line and dashed (red) line show the the A2796 and A2803 
Mgll absorption troughs. We find that the covering fraction largely 
determines the shape of the absorption trough. The troughs arc 
deeper at velocities where low-ionization gas covers a higher frac- 
tion of the continuum. 



3.2.1. Pure Partial-Covering Model 

An approach used extensively for quasar winds (Arav 
et al. 2001; Arav et al. 2005; Arav et al. 2008) is to assume 
that a fraction 1 — Cf(v) of the area of the continuum 
source is free of absorption at velocity v; and the optical 
depth distribution in front of the covered part of the 
source is constant. For this parameterization, the line 
intensities are given by 

I R (v)/I = l-C f (v) + C f e-^ (5) 



I B {v)/h = 1-Cf (v) + C f e- 2r(v \ (6) 

where t(v) is the optical depth at velocity v in the weaker 
line, here Mgn A2803. The coefficient in the argument of 
the exponential in Eqn.[S] depends on the relative oscilla- 
tor strengths and wavelengths of the doublet transitions. 
We binned the blue and red profiles to the spectral res- 
olution; and then computed the covering fraction 



C(v) 



I 2 

1 R 



2/ fl + l 



21, 



1 



(7) 



where the continuum has been normalized to unity, Iq — 
1. At velocities where measurement errors produced un- 
physical line itensities, Ir < Ir or I'i > Ir, we set the 



covering fraction to C(v) = 1 — Ir(v) or C(v) 
respectively. By substitution, the optical depth is 



t(v) = In 



C(v) 



I R {v) + C(v) - 1 



1.0, 



(8) 



Where statistical fluctuations result in Ir < Ib, mea- 
surement of the optical depth is not constrained. At 
velocities where Ir ~ Ir, both transitions must be op- 
tically thick; and this method cannot discriminate an 
optical depth of 3 from 30 in galaxy spectra of typical 
quality. The limited spectral resolution and low SNR of 
these data leave open the possibility of substantial sys- 
tematic errors in our best estimate t(v). We will refer 
to this method as the pure partial covering scenario with 
velocity-dependent covering fraction. 

Figure [6] shows the covering fraction of low-ionization 
absorption derived this way varies significantly across the 
trough. Our analysis necessarily neglects components 
narrower than 100 km s _1 , which may contribute to 
the total gas column but are unconstrained by our spec- 
tra. The covering fraction is near unity only where the 
trough is black. In the FSC0039-13 spectrum, this maxi- 
mum coverage is near the systemic velocity regardless of 
whether we remove a symmetric, zero- velocity absorption 
component before computing C(v). In FSC1630+15, the 
velocity of the low-ionization gas with the highest cov- 
ering fraction is between 300 and 400 km s _1 . The 
covering fraction falls steadily with increasing outflow 
speed in FSC0039-13, reaching 35% at 500 km s" 1 . 
Similarly, away from the velocity where the absorption 
trough presents minimum intensity, Cf(v) declines in 
FSC1630+15. The residual intensity correlates with the 
derived covering fraction. We conclude that the shape 
of the absorption trough, I(v), is strongly influened by 
the velocity dependence of the covering fraction of low- 
ionization gas in the outflow. 

Figure [7] compares the Mgn A2803 optical depth, t(v), 
to Cf(v). The low-ionization gas at velocity v ~ 
—600 km s _1 is apparently optically thick in the Mgll 
transitions. At higher speeds, the blending of the blue 
wing of the Mgn 2803 trough with the Mgn 2796 line 
prevents a unique comparison of Ir(v) to Ir(v). The 
solution to Eqn.[8]for for FSC0039-13 and FSC1630+15 
indicates that t(v) increases by a factor of ~ 2 and 1.5, 
respectively, as the covering fraction grows by a similar 
factor. Such a trend is physically plausible, even proba- 
ble, but is not required by the data. The method only 
returns a lower limit on optical depth, r>3. In our spec- 
tra, optical depths of 3, 30, and 300 yield essentially 
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identical line profiles; in contrast, damping wings (not 
seen) would distinguish r>10 4 . The lower limits on the 
optical depth, t(v), determines the lower limit on the 
ionic column density in the ouflow at a given velocity. 
Applying Eqn. 1 from Arav etal. (2001) to the optical 
depth in the Mgn A2803 transition, 

2l 3.7679 x 10 14 f , , , 

N[cm]= MMf l T{v)dv > (9) 

we estimate total ionic column densities of 7V(Mgn ) > 
6.4 x 10 14 cm" 2 and > 4.28 x 10 14 cm" 2 toward FSC0039- 
13 and FSC1630+15, respectively. 
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Fig. 7. — Optical depth of the Mgn A2803 line vs. Mgn covering 
fraction. We binned the Mgn absorption troughs to the spec- 
tral resolution and then computed Cf(v) and t(v) from Equa- 
tions[7]and [8] Triangles and pentagons represent FSC0039-13 and 
FSC1630+15, respectively. Where statistical errors led to unphysi- 
cal intensity ratios, Ig > Ir, no upper bound on the optical depth 
could be obtained; and no results are shown. Removal of a symmet- 
ric, zero-velocity component in FSC0039-13 assigns zero covering 
fraction and optical depth to the outflow near the systemic veloc- 
ity; but the fit to the outflow portion of the trough barely changes 
(open symbols). The data allow, but do not require, regions of 
the absorption trough with higher covering fraction to have higher 
optical depth. 



3.2.2. Inhomogeneous Absorption 

The pure partial covering scenario with velocity- 
dependent covering fraction is just one particular two- 
parameter (Cf, t) description of the relative intensities 
of the doublet troughs. With just two absorption troughs 
to fit, other two-parameter models would also exactly fit 
the data. A real test of the model requires spectra cov- 
ering additional transitions from the same ion, which is 
not possible with our data. We can, however, calculate 
parameter values for other models and examine how they 
change our view of the distribution of low-ionization gas 
in the outflow. 

One such model is inhomogeneous absorption across 
the source. This is appealing for AGN outflows because 
sharp edges in the spatial distribution of absorbing gas 
appear unrealistic because the gas lies at distances thou- 
sands of times greater than the size of the AGN emis- 
sion source (de Kool etal. 2002). Arav etal. (2005) 
present a simple model, where the optical depth distribu- 
tion across a uniform surface brightness source varies as 



t{v,x) = To(v)x a , where x describes the location (in the 
plane of the sky) across the projected area of the source; 
and x is confined to the range [0,1]. To produce similar 
intensities in the Mgn A2796 and A2803 lines, they re- 
quire the function t(x) be steep, e.g. a>8. A high value 
of the exponent works because it yields a factor of two 
change in the optical depth over a relatively small num- 
ber of spatial elements, x, resulting in I(2t) w I(t). For 
even moderately large values of t , the optical depth will 
be low over at least a few tenths of the source area, and 
the absorption troughs will not be black. 

In the ULIRG outflows, Ir ~ Ib over a wide veloc- 
ity range; but the Mgn troughs are nearly black over 
a smaller velocity range. The inhomogeneous absorber 
model requires unphysically large optical depths to repro- 
duce these Mgn absorption troughs. Starburst outflows 
also do not present the problem with pure-partial cover- 
ing as do AGN outflows. The rotation observed across 
some low-ionization outflows suggest the outflowing gas 
lies within a few kpc of the galaxy (Martin 2006). Hence, 
we find no reason to adopt an inhomogeneous-absorption 
model for ULRIG outflows. 

3.3. Fitting Velocity Components to the Absorption 
Troughs 

The direct method of the previous section cannot be 
applied to the absorption troughs in which the doublet 
transitions blend together and/or the absorption trough 
blends with one or more emission lines. A quantitative 
description of these absorption troughs requires fitting 
the parameters of functions describing each spectral com- 
ponent. 

A covering fraction, Doppler parameter, optical depth 
at line center, and Doppler shift describe the attenua- 
tion of the continuum by one component, or cloud, at 
every wavlength, Ij(X)/I = 1 - Cj + C^e" 7 ^. When 
the absorption trough is a blended doublet, the optical 
depth at a given wavelength may include significant con- 
tributions from each transition, and we add the doublet 
optical depths at that wavelength. At any particular ve- 
locity, the two transitions of the doublet must have the 
same C/, b, and Doppler shift; and atomic physics fixes 
the ratio of their optical depths. At the systemic veloc- 
ity of the Hel 5877 and Mgn emission lines, we fitted 
the amplitude of Gaussian intensity profile with velocity 
width matched to that measured for Ha . 

The relative intensities of the Mgn and Fen doublet 
troughs constrain these transitions to be optically thick, 
setting a firm lower limit on the optical depth, tq>3, at 
line center. The blue skew of the absorption troughs re- 
quires at least two velocity components. When we fitted 
such models, however, we found a number of parameter 
degeneracies. We could describe the velocity width of 
each trough by adding additional velocity components, 
increasing the Doppler parameter, or further increasing 
the optical depth of individual components. For exam- 
ple, ro can be made arbitrarily large by simply allowing 
the Doppler parameter b to shrink; and the product 6to of 
a component can be made smaller (or larger) by adding 
(or removing) additional velocity components. Previous 
studies (Rupke etal. 2005a; Martin 2005, 2006) chose 
the minimum number of velocity components required 
to describe the Nai absorption troughs. 

Simulations of winds show that a single sightline inter- 
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Fig. 8. — Fitted covering fraction of absorbing gas in each velocity component (red) and, when required, emission lines (green). Left 
Column: Fitted Mgn absorption troughs. All velocity components have the Doppler parameter fixed at b = 50 km s — 1 and an optical 
depth at line center of to. Middle Column: Velocity components fitted to the Mgl absorption troughs. The Doppler shift and b of each 
component matches the corresponding Mgn component, and we fitted the velocity-dependent covering fraction Cf(v) and minimum To- 
Right Column: Velocity components fitted to blended Nal doublet trough. The Doppler shift and b value of each component matches the 
corresponding Mgn component; we fit a velocity-dependent covering fraction Cf(v) and the minimum to. Dotted line shows the relative 
intensity of sky emission; an additive constant of 0.4 has been applied to the normalized sky intensity for display purposes. Although the 
fitted models are not unique; they demonstrate that all three line profiles can be described with the same kinematic components. 
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sects multiple shell fragments (Fujita etal. 2009). And 
high resolution spectra of starburst galaxies resolve some 
broad absorption troughs into components (Schwartz & 
Martin 2004). With this picture in mind, we assume 
a sightline intersects a number of velocity components, 
where these structures partially overlap spatially. To al- 
low for the possibility that some of the sightlines within 
the beam intersect different structures, we adopt the par- 
tial overlap model introduced by Rupke etal. (2005a). 
Each of the n successive velocity components attenu- 
ates the continuum (and the preceeding components) 
such that I(v) = YYj=i Ij( v )- A Maxwellian distribu- 
tion provides a reasonable description of the atomic ve- 
locities within each velocity component, so a Gaussian 
function models the optical depth distribution, r(A) = 

Toe~ ( - A ~ Ao ' ) /( A ° b / C ) , of each component. The covering 
fraction is constant across a given component but varies 
from one component to the next. The Doppler parame- 
ter, b, describes the turbulent motion along the line-of- 
sight in a given gas cloud. For purposes of illustration, we 
fix the Doppler parameter at b = 50 km s _1 for all veloc- 
ity components. For a typical minimum value of To ~ 3, 
this choice yields a line width at half the trough depth 
around 100 km s _1 , consistent with the spectral resolu- 
tion. We add velocity components until the fit statistic 
stops improving or a component becomes weaker than 
the detection limit. With this Doppler parameter, we 
need 6 velocity components to fit the Mgn absorption 
trough in FSC 0039-13 and 5 velocity components for 
the other 4 galaxies. 

Allowing To to vary among velocity components re- 
vealed no obvious trends with outflow velocity. The op- 
tical depth at the center of each component must be 
greater than three for the doublet transitions, but the 
probability distribution is very asymmetric with a signif- 
icant tail reaching values ~ 50 times higher than the min- 
imum value. (See the Monte Carlo approach described 
by Sato et al. (2009) for describing uncertainties in op- 
tical depth.) We obtained an equally good fit statis- 
tic when we required a single value of r describe all 
the velocity components. This assumption of a velocity- 
independent optical depth is our most dubious prior, and 
we discuss how to improve upon it in Section 13.3.31 We 
adopt this approach for the purposes of comparing the 
velocity-dependence of the gas covering fraction among 
objects and among transitions. For optically-thick veloc- 
ity components, the exact value o/tq(v) has little effect on 
Cf(v). Our approach assumes the velocity- depdendence 
of the covering fraction determines the trough shape. 

We first fit Cf(v) and the minimum value of the opti- 
cal depth to the Mgil absorption troughs. We found that 
the same kinematic components described the Fell and 
Mgi absorption troughs well. In Section 13.3. 1[ we show 
that the minimum optical depth in the weak lines raises 
our best estimate of the minimum optical depth in Mgn . 
We test the hypothesis these velocity components also 
describe the blended, Nal absorption troughs and com- 
pare covering fractions among species. Table [3] summa- 
rizes the fitted parameters. Figure [8] shows the resulting 
fits for Mgn , Mgi , and Nal . We plot Fen separately 
in Figure [9] 

3.3.1. Fitted Optical Depth 



For Mgn , we can compare the ionic column inferred 
from the fitted partial-overlap model to those estimated 
in Section 13.2.11 With b = 50 km s _1 , we construct 
the same ionic columns by setting ro(MgII 2803) = 2.75, 
which implies an integrated optical depth across each 



velocity component of r 
implied column densities, 



v/tFc-Vq&Ao = 2.3A. The 



Mcm-1 = L13 X 



(10) 



add up to iV(MgIl ) = 6.4 x 10 14 cm" 2 for the 6 
components describing FSC 0039-13. This model has 
the minimum tq that comes close to describing the 
Mgn absorption trough given our prior on b. The re- 
duced chi-squared statistic improves significantly, from 
1.425 to 1.143, when to is increased to 5. The fit statis- 
tic does not get much better as To is increased up to 
~ 50, so we consider to ~ 5 to be our lower limit when 
the Mgn absorption is considered independently. The 
pure partial-covering method also provided only a lower 
bound on To(Mgn ) due to saturation and low spectral 
SNR. 

Based on inspection of the Fen absorption troughs, we 
argue that the actual Mgn 2803 optical depths are 2 to 
3 times larger than the minimum required to fit the dou- 
blet. For a cosmic abundance ratio of gas-phase iron to 
magnesium, the optical depth in the stronger iron line 
is T (Feii 2600) - 0.57To(Mgn 2803)x(Feii )/x(Mgli ). 5 
For the weaker Fell line, To (Fen 2587) will always 
be 3.5 times lower than To(Feii 2600). It follows 
that To (Fen 2587) will be slightly less than unity if 
To(Mgn 2803) = 5 and if our estimate for the relative 
ionization correction, x(F eI1 ) ~ x(MgII ), is correct. 
This model predicts the Fen (2587) absorption trough 
is much shallower than the Fen (2600) trough as illus- 
trated for FSC 0039-13 and FSC 1630+15 in panels a 
and e, respectively, of Figure [HI 

For FSC 0039-13, this model spectrum is inconsis- 
tent with the shallow Fen (2587) trough, which requires 
To (Fen 2587) ~ 3 or larger. As addressed further in Sec- 
tion 13.41 a difference this large in ionization fraction is 
unlikely given predictions from photoionization model- 
ing. We conclude that the actual central optical depth 
in Mgn 2803 is closer to 18 (than 5). Panel d of Figured] 
shows that the resulting Mgn fit is really indistinguish- 
able from that presented in panel b. We find the fit to 
Fen 2587 in panel c preferable to that in panel a. 

For FSC 1630+15, the Fen 2587 absorption trough is 
slightly deeper than the prediction from the minimum 
T (Mgn 2803) model. Increasing T (Mgn 2803) and 
To (Fen 2587) to 10 and 1.6 provides an acceptable fit; 
larger optical depths make the Fen 2587 trough deeper 
than observed in this case. We examined the UVBLUE 
spectra of synthesized stellar populations to determine if 
other lines might be blended with the Fen 2600 line. We 
identified resonance transitions from Mnii , marked in 
Fig. [3] The Mnii 2594.499 transition does not cause the 

5 Although Fe and Mg have similar 1st and 2nd ionization poten- 
tials, the ionization equilibrium is complicated by the larger role of 
dielectronic recombinations for Fe (Shull & Van Steenberg 1982). 
The Murray et al. (2007) photoionization calculation for outflows 
indicates that the largest difference in the Mgn and Fell ionization 
fractions occurs at high gas density, 6 X 10 5 cm -3 , for the softest 
SEDs, when x(Fell ) can be twice as large as x(Mgll ). 
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Fig. 9. — Constraints on optical depth from comparison of Fen and Mgn absorption troughs. The Fell model in the left panel has 
t (Fen 2600) = 0.566T (Mgn 2803). The top and bottom two panels show FSC 0039-13 and FSC 1407+05, respectively. (a,b) FSC 0039- 
13. Assumes the minimum Mgn optical depth that fits Mgn 2795 and 2803, i.e. T (Mgn 2803) = 5. With to (Fen 2600) = 3, the Fen 2587 
lines are not nearly strong enough. (c,d) FSC 0039-13. Fitting the Fell 2587 trough requires To(Fen 2600) of at least 10, even when a 7th 
velocity component is included. Increasing To(Mgll 2803) to 17.6 is not inconsistent with the Mgn troughs. (e,f) FSC 1407+05. Assumes 
the minimum Mgn optical depth that fits Mgn 2795 and 2803, i.e. T (Mgn 2803) = 5. With T (Fell 2600) = 3, the Fen 2587 lines are 
not quite strong enough. The green line shows the fitted emission component. (g,h) FSC 1407+05. Fitting the Fell 2587 trough requires 
to (Fen 2600) of at least 5.7, even when a 7th velocity component is included. Increasing to (Mgn 2803) to 10 is not inconsistent with the 
Mgn troughs. These results suggest the Mgn optical depth is likely a factor of ~ 2 to 3 larger than the minimum required by the data. 
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depression seen at slightly shorter wavelengths because 
the stronger Mnii 2576.877 line is not detected. The 
SNR around the Fen 2587, 2600 lines is not adequate 
for further comparison in the other objects. 

Our observation of a single Mgi transition does not di- 
rectly constrain the Mgi 2853 optical depth. The optical 
depth in Mgi is 

r (Mgi 2853) = 6.1r (Mgn 2803)x(Mgi )/x(Mgn $11) 

We argued that r (Mgn ) reaches at least ~ 10 — 18, so 
the neutral Mg fraction must be less than 1 to 2 per- 
cent to produce an optically thin Mgi 2853 trough. It 
seems more likely to us that Mgi 2853 is optically thick. 
The covering fraction then determines the shape of the 
Mgi absorption trough, I(v) w 1 — Cf(v). Since cover- 
ing fraction dictates both the Mgi and Mgn absorption 
trough shape, they have a similar appearance, as shown 
quantitatively in the middle column of Figure [H That 
the Mgn troughs are deeper indicates that only a frac- 
tion of the cloud volume contains much neutral Mg. 

The fitted models in Figure [5] assume equiva- 
lent kinematic components for Nai and Mgn . 
For a cosmic abundance ratio, N(Mg)/N(Na) — 
19.07, the Nai optical depth, r (Nai 5898) = 
0.115ro(Mgn 2803)x(Nai )/x(Mgll ), will be lower 
than that of Mgn 2803 for similar ionization fraction. In 
Section [331 we argue that x(Nal )/x(Mgn ) may be of 
order unity for soft spectral energy distributions (SED) 
with L^ v I L^ IR ~ 10 -9 but will vary among outflows, 
dropping to x(Nai ) = 0.001 for hard starburst SEDs. 

When the weaker Nai line is optically thin, the fit- 
ted Nai models fall well short of the observed intensity 
in the redder half of the absorption trough. The mod- 
els shown in Figure [5] have To (5898) = 3 for 4 galax- 
ies and 5 for FSC 2349+24. The ratio of these lower 
limits to those for ro(Mgll 2803) in the same objects is 
0.3 or more. Since the photoionization models, and con- 
sideration of depletion factors, disfavor a situation with 
x(Nai )/x(Mgn ) > 1, the moderate optical thickness 
required for Nai suggests To(MgII 2803) could easily be 
a factor ~ 2 higher than that required by To (Fen 2587). 
As discussed in Section 3.1.1, however, the Nai and 
Mgn absorption troughs could probe physically distinct 
regions of the outflow. 

3.3.2. Constraints on Covering Fraction 

The attraction of forcing the same kinematic compo- 
nents is that we can directly compare covering fraction of 
different species as a function of velocity, as summarized 
in Tabled The fitted Cf values for the Fen velocity com- 
ponents come out similar to the values for Mgn . The 
maximum covering fraction in low-ionization gas occurs 
at minimum intensity by construction, and we again as- 
sociate the velocity of maximum covering fraction with 
the speed of swept-up shells at breakout. The covering 
fraction decreases towards higher outflow velocity. 

Independent of any particular model, the shallow 
troughs of Mgi relative to Mgn require a lower cover- 
ing factor for the former at every velocity. The xt~ 
fitting approach demonstrates that the shape of the 
Mgi absorption trough can be well described by the same 
velocity components as the Mgn troughs and further 
quantifies the change in covering fraction. Figure [TO] 
compares the fitted covering fractions for Mgn and Mgi . 



The covering fraction for neutral gas remains roughly half 
that of the Mgn at all velocities. Given the similar kine- 
matics but constant offset in projected area, we suggest 
that the Mgi absorption originates in denser regions of 
larger structures traced by Mgn . 

The fitted Nai covering fraction is similar to that fitted 
to Mgi for each velocity component. Sodium and mag- 
nesium have similar first ionization potentials, 5.1 eV for 
Na, and 7.6 eV for Mg. In the FSC0039-13 spectrum, 
emission filling could explain shallow Nai 5892 trough, 
relative to Mgi ; but the emission is not broad enough 
to account for the shallow Nai 5898 trough. We sug- 
gested that the FSC 0039-13 sightline intersects more 
of the galactic, gas disk than our other observations; 
and the high inclination of this disk means we probe the 
swept-up shell where it has stalled in the disk. The ex- 
tra attenuation at high inclination would likely produce 
larger than average differences between near-UV and op- 
tical continuum morhpology. High-resolution imaging 
in these bands may yield further insight into why the 
covering fraction of neutral, alkalai metals differs by al- 
most a factor of two between the near-UV and optical 
sightlines. In more typical ULIRGs, the similar kine- 
matics of Nai absorption and that of Mgi , suggest the 
Mgi velocity components and Cf(v) might be applied as 
template for deblending the Nai absorption trough into 
its Nai 5892, 98 components. 

3.3.3. Minimum Mass-Loss Rate 
In a spherical outflow geometry, the mass flux in the 
low-ionization outflow, M C (R) = flC f (R)v(R)p c (R)R 2 . 
The mass column, rhNuip) — J R1 p( r )dr, measured 
for any velocity carries a mass flux of M(vi) — 
rtrhNH(vi)Cf(vi)viRiR2/(R2 — Ri), where the distance 
and radial depth of the component are model dependent. 
We expect most of the mass-loading in the starburst re- 
gion and therefore that M(r) be approximately constant 
throughout the flow. For this to happen, the product 
T~Q,iCf t iVi needs to decrease as the radial distance to com- 
ponent i, r" 1 , increases. Our fitted values of covering 
fraction decline with increasing velocity, but the prod- 
uct ViCf(vi) grows with increasing outflow velocity. Ob- 
taining a constant mass flux requires the central optical 
depth to decline with increasing velocity. If we choose 
the minimum optical depth allowed for the highest veloc- 
ity component, then one can construct the values tq(v) 
fro v < max v needed to raise the mass flux at lower 
velocities to the same outflow rate, i.e. M(v) equal to a 
constant. We obtain the equivalent answer by estimating 
the mass flux from the highest velocity component. 

For FSC 0039-13, the velocity component 
at 571 km s _1 carries a mass flux of M > 
6 M yr- 1 x(Mgn )- 1 i?/Ai?(i?/10 kpc)(fi/47r), where 
the lower limit comes directly from the requirement that 
To(2587) > 2.7. The velocity component at 615 km s _1 in 
the Fell absorption troughs of FSC 1407+05 requires 
M w 5 M yr- 1 x(Feii )- l R/AR{R/W kpc)(fi/47f). In 
Section [33 we will argue that the ionization corrections 
for Mgn and Fen are of order unity. The unknown 
distance to the absorbing gas and radial thickness of the 
component leave order-of-magnitude uncertainty in the 
lower limit. 

Comparison of cosmic abundances and atomic data 
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Fig. 10. — Fitted covering fractions to Mgll (solid symbols) and Mgl (open symbols) absorption troughs in five ULIRGs. Each compo- 
nent's velocity is shown relative to the velocity of the component having minimum Mgll intensity; this component may be associated with 
the outflow velocity of a recently disrupted supershell. The shape of the absorption trough is dictated by the velocity-depedence of the 
covering fraction because each velocity component is optically thick. This modeling approach quantifies both the lower covering fraction 
of neutral Mg relative to singly-ionized Mg in the outflow and the trend towards lower covering fractions at high velocity. 
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indicates tq = 3 will yield the highest mass-loss 
rate for Nai 5898, with Fen 2587 about 66% as 
large. To make the implied mass columns agree, the 
minimum optical depth in Fen 2587 would need to 
be 4.5C/(Nai )/C/(Fell )x(FelI )/x(Nal ), or roughly 
2.25x(Feii )/x(Nai ). In Section we find the upper 
bound on the correction to the minimum optical depth, 
l/x(Nal ) could reach 10 3 for those ULIRGs with rela- 
tively hard SEDs. The lower limits on the mass fluxes 
likely yield an accurate picture for some systems but un- 
derestimate the true mass-flux in low-ionization gas by 
orders of magnitude in others. 

3.4. Gas Volume Density 

The volume density of the low-ionization gas varies 
widely among outflow models. Figures 4 and 5 of Mur- 
ray et al. (2007) illustrate the sensitivity of the Mgn and 
Nai ionization fractions to the gas density. Since the sat- 
urated lines provide only lower limits on the column den- 
sities of Mgi and Mgn , we cannot directly estimate the 
ionization fraction, thereby constraining the volume den- 
sity. The null-detection of absorption from collisionally- 
excited levels above the ground state provides a useful 
upper limit on the electron density, however. The most 
relevant excited transitions of Fen* are marked in Fig- 
ure 02 see Korista etal. (2008) for a summary. We com- 
pute this upper limit and then discuss how it constrains 
ionization fraction in the outflow. 

Above the critical density, the level populations ap- 
proach their Boltzmann ratio, and we can easily calculate 
the relative strengths of absorption lines from different 
energy levels. The lowest energy level above the ground 
state, E = 385 cm -1 , has the lowest critical density. Our 
spectra cover the second strongest line from the multiplct 
with this lower energy level, Fell* 2612.6542. 

In the FSC0039-13 spectrum, we expect 
Fen* 2612.6542 to be the strongest excited line. 
We place an upper limit on the observed equivalent 
width using the continuum SNR and assuming a 
line-width comparable to a spectral resolution cl- 
ement. The 5cr upper limit on the corresponding 
rest-frame equivalent width is 0.40A. For comparison, 
we also fitted a single-absorption-line model with 
b = 50 km s _1 and Ct = 1 to the spectrum. At the 
systemic velocity, we find r (Feil* 2612) < 0.30, a 
value typical of this spectral bandpass. We obtain a 
more conservative limit if we repeat this procedure 
where the largest continuum deviation occurs at - 
850 km s" 1 ; and we find r (Feii* 2612) < 0.62. This 
more conservative limit corresponds to a column density, 
iV(Feii* ) < 6.26 x 10 13 cm~ 2 per component. 

The optical depth in the individual velocity com- 
ponents of the weakest, detected resonance line is 
To(Feii 2587) > 3. Each of the 6 velocity components 
had b = 50 km s _1 , so the lower limit on the column den- 
sity in the ground state is 7V(Feii ) > 5.61 x 10 14 cm -2 . 
Our measurements place a solid upper limit on the rel- 
ative column densities in the first-excited and ground 
states, iV(Fell* )/iV(Feii ) < 0.11, and likely less than 
0.05. For densities well above the critical density this 
ratio would be ~ 0.75 with a slight dependence on tem- 
perature. Using the level population calculations from 
Figure 3 of Korista etal. (2008), we conservatively limit 
the electron density to logn e < 3.5 (or 3.4) for a tem- 



perature of 1 — 1.5 x 10 4 K (or 5 x 10 3 K), respectively. 
The stronger limit of N(Fcu* )/iV(Feii ) < 0.05, which 
applies if the highest density gas is at low velocity, lowers 
logn e to 3.1 (or 3.0), respectively, at T = 1-1.5 x 10 4 K 
(or T = 500 K). 

The gas density is important for understanding the re- 
lationship between the hot wind and the outflow observed 
in UV-optical absorption lines. Photoionization equilib- 
rium likely sets the temperature of the low-ionization gas 
at about 10 4 K. The density of the cool, low-ionization 
gas must be > 50 cm~ 3 to be in pressure equilibrium 
with the hot wind, where P^>(10 7 K)(0.05 cm" 3 ) ~ 
5 x 10 5 K cm~ 3 . Outflows accelerated by radiation pres- 
sure on dust grains do not require a hot wind at all (Mur- 
ray etal. 2005), allowing the low-ionization gas to have 
much lower density. Our upper limit on the gas density 
does not challenge the multi-phase models. 

Our result does eliminate the Murray et al. (2007) ion- 
ization models with n ~ 6 x 10 4 cm~ 3 or larger (see 
their Figures 4 and 5). We can confidently claim that 
Mgn is the dominant ionization state of Mg. For any 
reasonable starburst spectral energy distribution (SED), 
x(Mgn ) > 0.7; and the neutral fraction is less than 30%. 
The correction from the ionic columns of 7V(Mgn ) (and 
AT (Fen )) to the elemental columns will be relatively mi- 
nor. This knowledge of the Mg ionization balance sug- 
gests the Mgi 2853 optical depth significantly exceeds its 
minimum value of 3. Obtaining the same total column 
of Mg from the Mgi 2853 and Mgn 2803 components re- 
quires 



n- OSWK\ fi 1 x(Mgi) C / (Mgn) _ 
ro(2853) = 6.1 x(Mgn)C/(Mgi) r (2803), 



(12) 



or r (2853) up to 5.2r (2803) when x(MgI ) = 0.3 and 
C/(Mgl ) = 0.5C/(Mgn ). 

Eliminating very high density rules out a high neutral 
fraction of Na. Interpolating between the curves in Fig- 
ure 4 of Murray etal. (2007) to a density of 1000 cm~ 3 , 
the Nai ionization fraction plumments from unity, when 
L UV /L IR < IQ- 9 , to x(Nai ) - 10~ 3 when the SED 
hardens to Luv/LiR ~ 10~ 4 . In contrast, only at very 
high gas density does the neutral Mg fraction change 
rapidly with spectral hardness. The ionization potential 
of Mg is 2.5 eV higher than that of Na.. 

3.5. Emission 

Our spectrum of FSC2349+24 very clearly shows 
Mgn emission. The emission line appears to be a bit 
rcdshifted, but we argue that the Mgn absorption atten- 
uates the blue side of the emission profile. The only other 
spectrum that definitely presents Mgn is FSC1407+05. 
The Mgn 2803 profile shows the emission from to 
+250 km s _1 . The corresponding Mgn 2796 line is 
less obvious because of Mgn 2803 absorption at the 
same wavelength. We suggested that emission filling of 
the Mgn trough explains the deeper absorption in the 
Fell trough relative to the Mgn trough near systemic ve- 
locity. 

Weiner et al. (2009) discovered Mgn emission in a sub- 
set of z ~ 1.4 galaxies with outflows and argued that 
the composite spectrum of the non-emission galaxies pre- 
sented weak emission. They demonstated the presence 
of the emission by fitting the red-shifted portion of the 
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Mgn 2803 absorption trough with a symmetric absorp- 
tion component at the systemic velocity. Removal of this 
component yielded a redshifted emission component in 
Mgn 2796. This technique cannot be applied to four of 
five galaxies in our sample because their Mgn absorption 
troughs present no absorption at v > km s _1 . Remov- 
ing a symmetric, zero-velocity component fitted to the 
redshifted absorption trough in FSC0039-13 reveals no 
significant emission excess in Mgn 2796. 

The origin of the Mgn emission is not completely clear, 
but we adopt the hypothesis that the Mgn lines are ex- 
cited by recombination. The two objects presenting emis- 
sion are classified as Sey 2 on the basis of their optical 
emission-line ratios, suggesting the ionizing spectrum is 
harder. The absence of extended Mgn emission in the 
2D spectra of FSC1407+05 and FSC2349+24 rules out 
a scattering origin from a galaxy-scale nebula. 

The fitted profiles are shown in Figure [8] In 
FSC1009+47, the red side of the Fen troughs are 
marginally lower than those of Mgn , hinting at a hid- 
den emission component in the latter. However, inclu- 
sion of an emission component in the model for any of 
the other three galaxies does not improve the fit to the 
Mgn absorption trough. In fact, any emission compo- 
nent at the systemic velocity significantly degrades the 
fit if the maximum intensity is more than 10% of the 
continuum level. We conclude that only two of the 
five ULIRGS present significant Mgil emission. In con- 
trast, an He I 5876 emission improves the model for the 
Nal absorption trough in all five galaxies. The line is 
strongest in FSC0039-13, which is the only object classi- 
fied as an HII galaxy in our sample. 

3.6. Dynamical Ages of Targets 

In the popular "cool ULIRGs — ► warm ULIRGs — ► 
quasars" evolutionary scenario, as suggested by Sanders 
et al. (1988), one expects the AGN to provide the increase 
in bolometric luminosity. Advanced mergers in the 1 
Jy sample tend to have higher luminosity (Veilleux et al. 
2002). While our sample does contain rather luminous 
ULIRGs, their 25-60 fim colors, see Table [IJ are all red 
enough to be classified as cool ULIRGS, i.e. /25//6O < 
0.2. 

Following the morphological classification scheme out- 
lined in Veilleux et al. (2002), the most advanced mergers 
we observed are FSC0039-13 and FSC1407+05. These 
are single nuclei systems with a compact morphology 
and little tidal structure. The dynamically younger 
merger, FSC1009+47, presents tidal tails but the nu- 
clei have coalesced. The tails extend to 29 kpc in 
FSC1009+47 (Veilleux etal. 2002). Our second spec- 
trum for FSC1009+47 is extracted 44 kpc south (along 
the slit) of the nucleus. Veilleux et al. (2002) classify this 
feature as a prominent knot in a tidal arm rather than 
a second nucleus because it is not detected in their K' 
image. Pre-merger sources with separated nuclei are rare 
in the 1 Jy sample, but we observed at least two such ob- 
jects. FSC1630+15 is a close binary, separation 4.4 kpc 
(Veilleux et al. 2006); and FSC2349+24 is a wide binary, 
separation 14.5 kpc. 

Among this small sample we found no evolutionary 
trends in low-ionization outflow properties. We com- 
pared the outflow properties along the sightlines to each 
nucleus in the double systems. The outflow was always 



seen in each spectrum. 

4. DISCUSSION 

We detected outflows in 5 ULIRGs in ground-state ab- 
sorption from Nal , Mgi , Mgn , and Fen and illustrated 
the remarkably similar shapes of the absorption troughs 
in all these transitions. We then demonstrated that the 
velocity-dependence of the gas covering fraction deter- 
mines the trough shape. Similarities in trough shape 
suggest these species reside in the same, low-ionization 
gas structures. Because Mgn (andFeii ) absorption cov- 
ers a higher fraction of the continuum source, at a given 
outflow velocity, than do neutral Mg and Na, we argue 
that the absorbing clouds or filaments are not homoge- 
neous. 

The Doppler shift of the absorbing gas does not iden- 
tify its position along the line-of-sight. In particular, 
no consensus had been reached previously as to whether 
the highest velocity material detected resides close to the 
starburst region or at much larger radii. A close connec- 
tion likely exists between this neutral-atomic gas, out- 
flowing dust, and large-scale molecular outflows. Images 
revealing the presence of dust and molecular gas in M82 
out to a few kpc (Walter etal. 2002; Hoopes etal. 2005; 
Veilleux et al. 2009b) provide strong evidence for the sur- 
vival (or continuous creation) of dense clouds despite ab- 
lation by the hot wind and evaporation, as described (for 
example) by Marcolini etal. (2005). 

Along the outflow axis of M82, it remains unclear 
whether the denser material shares the same kinematics 
as the Ha emission. The Ha Doppler shift does increase 
with increasing distance along the minor axis (Heckman 
etal. 1990; Shopbell & Bland-Hawthorne 1998; Martin 
1998), consistent with acceleration. However, we lack 
such well-resolved, position- velocity information for most 
nearby (and all high-redshift) galactic outflows. In this 
section, we describe what the velocity-dependence of the 
gas covering fraction may imply about the location of the 
low-ionization absorbing gas along the line-of-sight. We 
begin by considering some pedantic dynamical models, 
starting with the blowout of a superbubble as described 
by DeYoung & Heckman (1994). 

4.1. Development of a Galactic Wind 

The thermalized energy from supernova explosions 
drives a shock front through a galaxy, sweeping interstel- 
lar gas into a thin, radiating supershell (Tenorio-Tagle 
& Bodenheimer 1988; Shull 1993). While the growing 
shell remains smaller than the pressure scale height of 
the interstellar medium, it plows through an ambient 
medium of essentially uniform density; and the shell de- 
celerates with velocity falling as v oc t~ 2 ' 5 , or equiva- 
lently v oc i?~ 2 / 3 , for a continuous injection of mechan- 
ical energy (Weaver etal. 1977). The shell mass grows 
linearly with the bubble volume at this stage, and the 
mass column through the shell grows linearly with ra- 
dius 

mN(r) = l/Srpo, (13) 

where po is the average density of the ISM. During the 
supershell phase, the covering factor of the low-ionization 
gas will be unity provided the shell traps the ionization 
front. 
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Emission-line images of nearby, starburst galaxies show 
that supershells outgrow their host galaxies. We expect 
the shell to accelerate when the density gradient becomes 
steeper than p(r) oc r~ 2 (McKee & Ostriker 1988), a 
highly unstable situation in which a dense shell pushes 
on more rarefied gas. Numerical simulations show that 
the shell breaks up due to hydrodynamic instabilities at 
a few pressure scale heights (MacLow, McCray, &: Nor- 
man 1989). The evolution of the shell fragments has 
received relatively little attention, but they are clearly 
one source of low-ionization gas that will absorb contin- 
uum emission. Additional sources of low-ionization gas 
include hydrodynamic instabilities induced by shear at 
the disk - wind interface (Heckman et al. 2000) and pre- 
existing, interstellar clouds over run by the supershell 
(Cooper etal. 2008). 

A primary origin for the low-ionization outflows in shell 
fragments is appealing because it naturally explains the 
velocity offset of the Mgn absorption troughs from the 
systemic velocity. Due to the high central concentration 
of gas in ULIRGs, the pressure scale-height is several 
times smaller than the value h z ~ 100 pc typical of nor- 
mal galaxies; and theblowout radius Rq ~ 3h z ~ 200 pc. 
By the time of blowout, the mean speed of a shell has 
dropped to 

/ T v 1/3 

v(R) = 164 km s" 1 ^ r 

V ' \7.08 x 10 43 erg s" 1 J 

where the mass per H atom is m = 1.4mjy in the ambi- 
ent medium and L w , the rate of mechanical energy injec- 
tion, has been scaled to a SFR representative of ULIRGs, 
i.e. 100 M Q yr -1 . 6 A shell velocity ~ 200 km s _1 can 
describe the Doppler shift of the deepest part of the 
Mgn absorption troughs in Figure lb-le. A sightline at 
high inclination (i.e. close to the plane of the gas disk) 
would intersect a lower velocity shell due to the higher 
average gas density in the plane. We appeal to this sce- 
nario as a plausible explanation for the kinematics of the 
Mgn absorption troughs in FSC0039-13, which are deep- 
est near the systemic velocity. 

Numerical simulations of ULIRGs (e.g. Fujita et al. 
2009) confirm that shell velocities decline to 200- 
300 km s _1 by the time of blowout and suggest that, 
following blowout, the hot interior of the bubble acceler- 
ates outward creating a hot wind with terminal velocity 

Vh ~ V^3c s ~ 940T 7 1 ^ 2 km s _1 . Fujita etal. did not 
follow the shell fragments very far into the halo; they 
simply assumed the clouds coast outwards on ballistic 
trajectories. Whether or not cosmic rays (Breitschwerdt 
2008; Everett etal. 2008; Socrates etal. 2008), radiation 
pressure, and/or the ram pressure of the hot wind further 
accelerates the low-ionization gas remains an important 
question. 

Murray et al. (2005) analytically modeled low- 
ionization outflows accelerated by the ram pressure of a 
hot wind and outflows accelerated by radiation pressure. 
It remains challenging to observationally distinguish not 

6 Value of L-w from SB99 continuous star formation model (Lei- 
therer etal. 1999) for 1 M yr _1 in 1 to 100 M stars. 



only these two types of momentum-driven outflows but 
also the energy-conserving, ballistic trajectories. The for- 
mer accelerate quickly and then coast. The latter coast 
for a large distance before gradually decelerating. Fac- 
tors favoring radiative-driving include an empirical corre- 
lation between outflow speed and escape velocity (Martin 
2005), the high dust content of ULIRGs, and the high 
luminosity of ULIRGs. Reasonable objections include 
the small number of dwarf galaxies used in the outflow 
speed correlation, the difficulty of coupling radiative mo- 
mentum to the gas in dust-poor, lower luminosity dwarf 
galaxies, and evidence for the presence of hot winds in 
ULIRGs (Sciortino, V. & Martin, C. L. in preparation). 

For purposes of illustration, we consider acceleration 
by the ram pressure of a hot wind following blowout. 
The shell radius, i?o, and velocity, vq, at blowout set 
the initial conditions at the start of the wind phase. In a 
spherical outflow geometry, the ram pressure drops as the 
inverse-square of the radial distance, so clouds accelerate 
over a relatively short spatial scale. The terminal velocity 
of any shell fragment depends on its column density, with 
the low columns reaching, at most, the hot wind veloc- 
ity, v c w Vh (Martin 2005). The absorption troughs pre- 
sented in this paper require the covering fraction of shell 
fragments to decrease with increasing velocity. This situ- 
ation can be achieved by tuning the distribution of cloud 
column densities. The lowest column density fragments, 
which reach the highest terminal velocities, would need to 
be relatively rare, covering less area than higher column 
density clouds. Alternatively, the geometrical dilution of 
clouds offers a simpler way to achieve the desired result. 
We demonstrate this largely geometrical effect in Sec- 
tion 14.21 Further numerical work, beyond the scope of 
this paper, will be required to address the broader prob- 
lem. 

4.2. Geometrical Interpretation of the 
Velocity- Dependent Covering Fraction 

To illustrate the importance of geometrical dilution, let 
the locus of fractures in a shell define individual clouds. 
Suppose, for simplicity, these clouds all have the same 
column density and intitial velocity, equal to that of the 
shell just before blowout. If the area of a cloud does 
not change as it is moves outwards (think of a brick-like 
cloud), then the covering fraction of clouds decreases as 
C f (R) = C f (R Q )(R /R) 2 . 

Most shell fragments are unlikely self-gravitating struc- 
tures, however; and we expect them to expand as they 
fly outwards due to the drop in ambient pressure. The 
sound crossing time in a cloud is likely short enough 
for the cloud to maintain pressure equilibrium with the 
hot wind, P c ~ Ph- The clouds expand adiabatically, 
so we can calculate their size at any outflow radius, R, 
and compare the cloud area to the that of a solid shell, 
thereby estimating the covering fraction, Cf(R). 

Mass conservation in a steady-state, constant velocity 
hot wind requires the density to fall as oc r~ 2 . For an 
isothermal hot wind, the resulting increase in the volume 
of low-ionization clouds is 

V C (R)= (A) 2/7 V c (i? ), (15) 

where "f c = 5/3 for a monatomic, ideal gas. For roughly 
spherical clouds, the increase in cloud volume is accom- 
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panied by an increase in cloud area, A c oc V c . Defining 
the covering fraction as 



C f (R) _ Ac(R) 47ri?2 



C f (Ro) 4irR 2 AciRo)' 
these relations yield a covering fraction, 



R 



C f (R) _ 



4/3 7 c-2 



(16) 



(17) 



that falls as R~ 12 . 

If the hot wind cools adiabatically, then Th oc i?~ 4 / 3 ; 
and the pressure falls off faster with radius. The lower 
pressure allows the clouds to expand faster than in the 
isothermal case. Repeating the steps in the previous 
paragraph, we find 



CAR) 

C f (Ro) 



R 
Ro 



47h/37c-2 



(18) 



Assuming 7 = 5/3 for both the hot wind and the clouds, 
we find 



C f (R) 
Cf(Ro) 



R_ 

Ro 



-2/3 



(19) 



The low-ionization clouds cannot expand quickly enough 
to keep up with the geometrical dilution inherent to 
spherical, outflow geometry. Their covering fraction 
must decrease with increasing outflow radius. 

4.3. Is Acceleration of the Low- Ionization Gas 
Required? 

The absorption trough measurements require lower C / 
at higher velocity. Geometrical dilution produces lower 
covering fractions at larger outflow radii. Combining 
these two results implies that the gas producing the high- 
est velocity absorption resides at the largest distance 
from the starburst. While we find this argument illu- 
minating, some of the underlying assumptions should be 
examined. 

First, the outflow geometry is uncertain. It should be 
roughly spherical on spatial scales larger than the galaxy, 
but the outflow geometry may be better described as 
cylindrical at blowout. We know the absorbing gas lies 
close to the gaseous disk in some ULIRG outflows be- 
cause the outflows present rotation in Nal (Martin 2006). 
In a version of the above model with cylindrical geom- 
etry, both the covering fraction and column density of 
'brick-like' clouds become independent of height above 
the disk. The decrease in covering fraction would require 
other physical effects such as cloud ionization, evapora- 
tion, or ablation to become important at the higher out- 
flow velocities. 

Second, blowout may yield a distribution of cloud 
column densities and velocities, something three- 
dimensional numerical simulations may soon be able to 
address. The fastest fragments - whether determined by 
blowout, acceleration by the hot wind, or acceleration by 
the starburst radiation - will reach the largest distances. 
Our models indicate the low-ionization cloud reach their 
maximum distance long after the starburst activity has 
ceased. It follows that if we observe such outflows during 
the starburst phase, the highest velocity absorption will 



come from material at the largest radii. Our underly- 
ing explanation for the shape of the absorption troughs 
would remain geometrical dilution; however, acceleration 
of the low-ionization outflow would not be required be- 
yond blowout. 

4.4. Measurement of Terminal Velocities and Estimates 
of Spatial Extent 

Comparison of the Mgn , Mgi , and Nal absorption 
troughs indicates the highest velocity gas sometimes 
escapes detection in Nal and Mgi . For example, 
in FSC0039-13, we detect Mgi and Nal absorption 
to 300 km s _1 ; but the Mgll absorption demon- 
strates that the outflow persists to 500 km s _1 . In 
FSC1407+05, we measure terminal velocities of 350, 600, 
and 750 km s _1 from Nal , Mgi , and Mgn , respectively. 
We find our spectral signal-to-noise ratio insufficient to 
detect the higher velocity components of these outflows 
in Nal and Mgi . 

The Mgn lines allow more robust terminal velocity 
measurements for a number of reasons. First, the 
Mgn 2796, 2803 transitions have a higher value of 
N{X)/N H f than docs Nai 5892, 5898. Second, expected 
ionization corrections favor singly-ionized Mg, Fe, and 
Na over their neutral species. Third, we find a larger 
cloud covering fraction for the singly-ionized lines than 
for neutral lines. 

We emphasize that the highest velocity absorption de- 
tected in Mgll appears to be limited by covering frac- 
tion. At continuum SNR ~ 10, the absorption trough 
blends with the continuum where C/(u)<0.1, regardless 
of column density. The gas at the largest radii may well 
go undetected by typical absorption-line measurements. 
For example, if we assume unity covering fraction at Ro, 
then gas beyond R <~ 3.2i?o cannot be detected if C/ 
falls as R~ 2 . In the more realistic case, Cf oc i?~ 2 / 3 , 
the covering fraction drops to 0.1 at 32i? - For a launch 
radius R w 200 pc, our measurement would detect gas 
out ot 6.4 kpc. This distance exceeds the extent of the 
Ha filaments in many nearby starburst galaxies but re- 
mains considerably smaller than the recently detected 
soft, X-ray halos in ULIRGs (Sciortino & Martin, in 
prep.). Spectroscopy of background light sources at pro- 
jected separations of just a few tens of kpc from starburst 
and post-starburst galaxies should be more effective than 
linc-of-sight studies for determining the spatial extent of 
the low-ionization outflow, but such studies need to ac- 
count for the geometrical dilution of the clouds. 

5. SUMMARY AND CONCLUSIONS 

We presented the first comparison of optical and near- 
ultraviolet absorption troughs in ULIRG spectra. We 
detected outflows in 5 ULIRGs in ground-state absorp- 
tion from Nai , Mgi , Mgn , and Fen . Previous obser- 
vations of ULIRG outflows have been limited to the Nai 
doublet, and blending of the Nal 5890, 5896 lines compli- 
cated measurement of the absorption trough shape. We 
summarize the primary, empirical results. 

• Comparison of the unblended, doublet components 
of Mgn and Fen distinguish the effects of opti- 
cal depth and covering fraction in determining the 
shape of the absorption troughs. The high optical 
depth in these transitions at all outflow velocities 
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places a lower limit on the column density. The 
non-zero intensity requires partial coverage of the 
continuum source by the low-ionization gas over a 
broad velocity range. 

• The covering fraction in all four ions decreases as 
the outflow velocity increases beyond the velocity 
of minimum intensity, or equivalcntly the velocity 
of maximum covering fraction. At a given velocity 
the Mgi covering fraction is roughly half that mea- 
sured in Mgn (and Fell ), and the covering frac- 
tion in Nai is less than or equal to that measured 
in Mgi . 

• Accounting for these differences in ionic cover- 
ing fraction, and taking spectral SNR into con- 
sideration, we detect Mgn , Fen , Mgi , and 
Nai absorption over the same velocity range. The 
decrease in covering fraction with increasing ve- 
locity suppresses the absorption signature of the 
highest velocity gas. The higher covering frac- 
tion of Mgn (and Fen ) yields detection to higher 
velocity than indicated by Mgi (or Nai ). Any 
comparison of terminal velocities measured from 
Nai and Mgn must take the bias introduced by the 
velocity-dependent covering fraction into consider- 
ation. Many measurements exist for Nai at z ~ 
0.6, whereas ground-based, optical spectrographs 
can measure outflows with the Mgn doublet over 
the broad redshift range from ~ 0.25 to z ~ 2.5. 
Caution should be exercised when examining evo- 
lution in outflow properties between these samples. 
For completeness, we point out that our spectra do 
not cover high-ionization transitions like Ovi 1032, 
1038, which reveal higher-velocity gas in some star- 
burst outflows (Grimes etal. 2009). 

• The absence of Fell* lines, principally the A2612 
transition, place an upper limit on the volume den- 
sity of n e < 10 3 - 5 cm~ 3 , and likely < 10 31 cm' 3 . 
When present in quasar outflows, low-ionization 
absorption troughs from excited or metastable 
states indicate much higher volume density n e — 
10 4 - 4 cm~ 3 with less than 20% scatter (Korista 
etal. 2008; Arav etal. 2008). In the lower density 
ULIRG outflows, Mgn will be the dominant ioniza- 
tion state in the low-ionization outflow over a broad 
range of spectral hardness. In contrast, at these 
densities the Nai ionization fraction remains very 
sensitive to spectral hardness. We suggest that the 
harder radiation field in dwarf galaxies likely ex- 
plains the lower fraction of dwarf starburst outflows 
detected in Nai relative to ULIRGs, which are al- 
most always detected in Nai (Martin 2005; Rupkc 
etal. 2005b). 

• We found Mgn emission in two Sey 2 ULIRGs and 
He I emission in all spectra. 

These results provide new insight into the relationship 
of the low-ionization outflow and the hot wind. We de- 
fer physically-motivated models of the absorption trough 
shape to another paper (Martin 2009, in preparation)) 
but summarize key aspects of the emerging, physical pic- 
ture here. 



• We associate the velocity of maximum covering 
fraction with that of a swept-up shell of inter- 
stellar gas at the time of blowout. Factors mo- 
tivating this interpretation include the Doppler 
shift of the trough minimum (0 to -400 km s _1 ), 
the large width of the absorption troughs (up to 
800 km s _1 ), the column density lower limits, and 
results from recent numerical simulations (Fujita 
etal. 2009). Regardless of their physical origin, 
however, the large velocity width of the absorption 
troughs require contributions from multiple struc- 
tures along the sightline. The relative shapes of the 
doublet troughs require these clouds or filaments to 
transition sharply (spatially) from opaque to opti- 
cally thin gas. The neutral alkalai metals reside in 
the same kinematic structures as the singly ionized 
metals but fill a smaller fraction of that volume. 

• The most significant result of our study may be 
the discovery of a velocity-dependent covering frac- 
tion in low-ionization outflows. The simplest inter- 
pretation is geometrical. The dilution associated 
with the spherical expansion of a population of ab- 
sorbers causes their covering fraction to decrease 
with increasing radius. We showed that the adi- 
abatic expansion of clouds in pressure equilibrium 
with the hot wind is not fast enough to offest this 
dilution. In the context of this physical scenario, 
our result implies that the high-velocity gas detected 
in the absorption trough is at larger radii than the 
lower velocity (and higher covering fraction) gas. 
This mapping between velocity and relative radius 
indicates allows requires acceleration of the low- 
ionization gas. Either shell blowout or subsequent 
momentum-driving by a hot wind and/or radiation 
pressure could cause this acceleration. 

• In an alternative scenario, the shell fragments into 
clouds with a wide range of column densities. Due 
to momentum conservation, the ram pressure of 
the hot wind (or radiation pressure) would accel- 
erate the lowest column density fragments to the 
highest terminal velocities. This relationship might 
even arise from blowout alone without requiring 
a momentum-driven wind phase. Either way, the 
new empirical constraint, requiring lower covering 
fraction at higher velocity, would imply the cov- 
ering fraction of low-ionization gas increases with 
increasing column density. We find the first, geo- 
metrical explanation more appealing due to its sim- 
plicity. 

• The lower limits on Nai and Fell column density 
provide the highest (i.e. strongest) lower limits 
on the mass outflow rate in low-ionization gas. 
Better limits can be obtained by either determing 
the Nai ionization fraction or observing bluer, 
Fen transitions with lower oscillator strengths. 

• Assuming that photoionization equilibrium likely 
sets the temperature of the low-ionization outflow, 
T c ~ 10 4 K, the upper limit on the volume den- 
sity allows the pressure of the cold outflow to be 
as high as P c /fc<10 7 K cm~ 3 . Tighter constraints 
on the volume density would be valuable, as very 
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low density would be incompatible with pressure 
equilibrium with a hot wind, P c ~ Ph, and favor a 
purely radiatively-driven outflow. 

Our conclusion about the acceleration of the low- 
ionization outflow, while illustrative, is subject to an as- 
sumption about outflow geometry and a bias towards 
a shell origin for the low-ionization gas. Our empiri- 
cal results, however, clearly challenge plausible dynam- 
ical models. Acceptable models need to reproduce the 
velocity-dependence of the low-ionization covering frac- 
tion, explain the lower covering fraction of neutral, al- 
kalai metals relative to low-ionization species, and be 
consistent with the observed constaints on both column 
density and volume density. With sufficient spectral sen- 
sitivity, all of these outflow properties can be measured 



over a very broad redshift range, so we expect any evolu- 
tionary effects with galaxy mass and/or cosmic time to 
eventually be measured. 
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V (km/s) V (km/s) 

(a) LRIS spectra of FSC0039-13 

Fig. 1. — Normalized intensity vs. velocity, where the velocity is relative to the systemic velocity determined from emission lines. 
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(b) LRIS spectra of FSC1009+47 



Fig. 1. — continued 
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(b) LRIS spectra of FSC1407+05 



Fig. 1. — continued 
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(b) LRIS spectra of FSC1630+15 



Fig. 1. — continued 
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Fig. 1. — continued 
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TABLE 1 



Object 


f25/f60 


f60/fl00 


Spectral Type 


log(L/ij/ L ) 


SFR-S 


SFR-C 












( M yr- 1 ) 


( M Q yr- 1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


FSC00397-1312 


0.180 


0.963 


HII (V99) 


12.81 


1110 


653 


FSC10091+4704 


0.068 


0.761 


LINER(V99,K98b) 


12.57 


641 


377 


FSC14070+0525 


0.131 


0.797 


Sey2 (K98b) 


12.66 


789 


464 


FSC16300+1558 


0.047 


0.744 


LINER (K98b) 


12.53 


585 


344 


FSC23498+2423 


0.118 


0.703 


Sey2 (V99) 


12.31 


352 


207 



Note. - (1) ULIRG name. (2) f25/f60 color computed 
from far-infrared fluxes in KS98. (3) f60/fl00 color computed 
from far- infrared fluxes in [KS98]. (4) Spectral type from Kim 
ctal. (1998b) - FSC1009+47, FSC1407+05, FSC1630+15, - and 
Veilleux ct al. (1999) - FSC0039-13, FSC1009+47, FSC2349+24. 
(5) LiR calculated using the prescription of Sanders & Mirabel 
(1996). For FSC objects, the luminosities from Kim etal. (2002) 
and Veilleux et al. (1999), which assume luminosity distances based 
on Ho = 75 km s -1 Mpc -1 and go = 0, were converted to the 
cosmology used throughout this paper. (6) For Salpcter IMF from 
0.1 to 100 M , the star formation rate is SFR = L /fl /5.8xl0 9 L Q , 
where Lju is the bolometric luminosity (Kcnnicutt 1989). (7) SFR 
for Chabrier IMF, i.e. K89 SFR by 1.7. 



TABLE 2 



Object 


Rcdshift 


PA 


Res. 


SNR 


Date Observed 


Notes 








(deg) 


(km/s) 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




FSC00397-1312 


0.26171 


-54.0 


160 


9.4 


2007 Nov. 1 


Advanced Merger 




FSC10091+4704 


0.24508 


168.0 


130 


5.0 


2004 Jan. 26, Mar. 16-17 


Diffuse Merger 




FSC14070+0525 


0.26602 


-25.0 


110 


8.2 


2004 Jan. 26, Mar. 16-17 


Advanced Merger 




FSC16300+1558 


0.24200 


-45.0 


155 


6.2 


2004 Mar. 16-17, 2007 Oct. 6 


Pre-Merger, Separation 


4.4 kpc 


FSC23498+2423 


0.21249 


-45.9 


160 


6.4 


2007 Nov. 1 


Pre-Merger, Separation 


14 kpc 



Note. — (1) Target. (2) Rcdshift measured from emission lines 
including Ha , [Nil] 6584,48, [SII] 6717,31, and [OI] 6300,64. The 
absolute velocities have not been corrected to the Local Standard 
of Rest. (3) Position angle of the longslit measured east of north. 

(4) Spectral resolution denned by the measured full width at half 
maximum intensity of the arc lamp lines. If frames from different 
runs are combined, then the poorer of the two resolutions is listed. 

(5) SNR per pixel is measured in a lOOA bandpass immediately 
blueward of Mgll. (6) Date(s) observed. (7) The Veilleux etal. 
(2002) merger classification is compared to the structure observed 
along the longslit. 
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TABLE 3 
Fitted Covering Fractions 
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83 


0.74 + 0.03 


17.7 


1.136 


0.41 ± 0.22 
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1.417 


0.17 + 0.02 
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2.335 


0.99 ± 0.04 


2.8 


2.452 


0039-13 


-2.6 


0.83 ± 0.03 


17.7 


1.136 


0.57 ± 0.38 
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1.417 


0.04 ± 0.02 
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2.335 


0.90 ± 0.05 


2.8 


2.452 


0039-13 
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0.47 + 0.06 


17.7 


1.136 


0.20 ± 0.08 
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1.417 


0.11 + 0.02 
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2.335 


0.72 ± 0.05 


2.8 


2.452 


0039-13 
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0.45 ± 0.05 


17.7 


1.136 


0.24 ± 0.07 
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1.417 


0.05 ± 0.02 
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2.335 


0.50 ± 0.06 


2.8 


2.452 


0039-13 
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0.33 ± 0.02 


17.7 


1.136 


0.19 + 0.04 


4 


1.417 


0.09 + 0.01 


4 


2.335 


0.31 ± 0.04 


2.8 


2.452 


0039-13 
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0.27 + 0.02 


17.7 


1.136 


0.21 ± 0.03 
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1.417 


0.07 + 0.01 
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2.335 


0.43 ± 0.04 


2.8 


2.452 
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0.66 + 0.05 
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0.955 


0.17 + 0.11 


3.4 


1.243 


0.16 + 0.05 
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0.372 








1009+47 


-361 


0.55 + 0.05 


5 


0.955 


0.11 ± 0.08 


3.4 


1.243 


0.08 ± 0.05 


3 


0.372 








1009+47 


-239 


0.51 + 0.05 


5 


0.955 


0.25 + 0.10 


3.4 


1.243 


0.17 + 0.04 


3 


0.372 








1009+47 


-507 


0.51 + 0.04 


5 


0.955 


0.07 ± 0.07 


3.4 


1.243 


0.02 ± 0.04 


3 


0.372 








1009+47 


-692 


0.38 + 0.04 


5 


0.955 


0.10 ± 0.06 


3.4 


1.243 


0.04 ± 0.06 


3 


0.372 








1407+05 


-89 


0.29 + 0.04 


10 


0.989 


0.16 ± 0.08 


6.8 


1.045 


0.11 + 0.02 


3 


0.702 


0.52 ± 0.07 


1.6 


1.232 


1407+05 


-187 


0.40 ± 0.07 


10 


0.989 


0.08 ± 0.07 


6.8 


1.045 


0.10 + 0.02 


3 


0.702 


0.14 + 0.11 


1.6 


1.232 


1407+05 


-265 


0.62 + 0.04 


10 


0.989 


0.28 ± 0.09 


6.8 


1.045 


0.17 + 0.02 


3 


0.702 


0.55 ± 0.08 


1.6 


1.232 


1407+05 


-434 


0.67 + 0.04 


10 


0.989 


0.37 ± 0.05 


6.8 


1.045 


0.20 + 0.01 


3 


0.702 


0.49 ± 0.07 


1.6 


1.232 


1407+05 


-615 


0.51 + 0.03 


10 


0.989 


0.15 ± 0.04 


6.8 


1.045 


0.14 + 0.02 


3 


0.702 


0.59 ± 0.07 


1.6 


1.232 


1630+15 


-114 


0.28 + 0.05 


5 


1.086 


0.10 ± 0.09 


3.4 


1.014 


0.13 + 0.03 
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0.491 








1630+15 


-252 


0.81 + 0.05 


5 


1.086 


0.24 + 0.12 


3.4 


1.014 


0.21 + 0.03 


3 


0.491 








1630+15 


-362 


0.82 + 0.05 


5 


1.086 


0.42 ± 0.14 


3.4 


1.014 


0.16 + 0.04 


3 


0.491 








1630+15 


-491 


0.60 ± 0.05 


5 


1.086 


0.22 ± 0.07 


3.4 


1.014 


0.12 + 0.03 


3 


0.491 








1630+15 


-676 


0.21 + 0.04 


5 


1.086 


0.08 ± 0.05 


3.4 


1.014 


0.07 + 0.04 


3 


0.491 








2349+24 


-71 


0.27 + 0.03 


2 


2.247 


0.12 ± 0.07 


1.4 


2.254 


0.08 + 0.007 


5 


4.259 








2349+24 


-182 


0.28 + 0.03 


2 


2.247 


0.04 ± 0.05 


1.4 


2.254 


0.05 ± 0.007 


5 


4.259 








2349+24 


-298 


0.38 ± 0.03 


2 


2.247 


0.30 ± 0.07 


1.4 


2.254 


0.02 ± 0.007 


5 


4.259 








2349+24 


-423 


0.27 + 0.03 


2 


2.247 


0.27 ± 0.06 


1.4 


2.254 


0.02 ± 0.007 


5 


4.259 








2349+24 


-591 


0.17 + 0.03 


2 


2.247 


0.17 + 0.05 


1.4 


2.254 


0.008 + 0.007 


5 


4.259 









NOTE. — (1) Galaxy. (2) Dopplcr shift of component in km s _1 . All components have a Dopplcr parameter b = 50 km s _1 , or 650 = 1. (3) Fitted covering fraction for Mgll (4) Minimum 
optical depth at line center for Mgll 2803. The lower limit on the ionic column density is 7V(Mgll ) > 3.90 X 10 13 cm _2 Tob5o- (5) Fit statistic for the Mgll doublet. (6) Fitted covering fraction for 
Mgl (7) Minimum optical depth at line center for Mgl 2853. The lower limit on the ionic column density is iV(Mgl ) > 6.39 X 10 12 cm - 2 to&5o- (8) Fit statistic for the Mgl line. (9) Fitted covering 
fraction for Nal . (10) Minimum optical depth at line center for Nal 5898. The lower limit on the ionic column density is 7V(Nal ) > 1.78 X 10 13 cm - 2 ro??5o- (H) Fit statistic for the Nal doublet. 
(12) Fitted covering fraction for Fell . (13) Minimum optical depth at line center for Fell 2600. The lower limit on the ionic column density is 7V(Fcll ) > 1.87 X 10 14 cm - 2 ro^5o- (14) Fit statistic 
for the Fell doublet. 



